Astrophysical Journal 



A Spitzer IRS Low-Resolution Spectroscopic Search for Buried 
AGNs in Nearby Ultraluminous Infrared Galaxies — A Constraint 
on Geometry between Energy Sources and Dust — 

O ■ Masatoshi Imanishi 

O 

^ ■ National Astronomical Observatory, 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 
^\ C. C. Dudley 



X 



Naval Research Laboratory, Remote Sensing Division, Code 7211, Building 2, Room 24OB, 
4555 Overlook Ave SW, Washington DC 20375-5351, U.S.A. 



(N 
> 

en . Roberto Maiolino 

^ ■ INAF - Osservatorio Astronomico di Roma, via di Frascati 33, I-OOO4O Monte Porzio 

^ ■ Catone, Roma, Italy 



Philip R. Maloney 



O \ formerly Center for Astrophysics and Space Astronomy, University of Colorado, Boulder, 

^ ■ CO 80309-0389, U.S.A. 

Takao Nakagawa 

Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 3-1-1 
Yoshinodai, Sagamihara, Kanagawa 229-8510, Japan 

and 

Guido Risaliti^ 

INAF - Osservatorio Astrofisico di Arcetri, L.go E. Fermi 5, 1-50125 Firenze, Italy 

ABSTRACT 



^Department of Astronomy, School of Science, Graduate University for Advanced Studies, Mitaka, Tokyo 
181-8588 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, U.S.A. 



- 2 - 



We present the results of Spitzer IRS low-resolution infrared 5-35 /im spec- 
troscopy of nearby ultraluminous infrared galaxies (ULIRGs) at z < 0.15. We 
focus on the search for the signatures of buried active galactic nuclei (AGNs) 
in the complete sample of ULIRGs classified optically as non-Seyferts (LINERs 
or Hll-regions). In addition to polycyclic aromatic hydrocarbon (PAH) emission 
features at 6.2 /im, 7.7 /im, and 11.3 fim, the conventional tool of starburst-AGN 
separation, we use the optical depths of the 9.7 /xm and 18 fim silicate dust ab- 
sorption features to infer the geometry of energy sources and dust at the nuclei of 
these ULIRGs, namely, whether the energy sources are spatially well mixed with 
dust (a normal starburst) or are more centrally concentrated than the dust (a 
buried AGN). Infrared spectra of at least 30%, and possibly 50%, of the observed 
optical non-Seyfert ULIRGs are naturally explained by emission consisting of (1) 
energetically insignificant, modestly obscured (Ay < 20-30 mag) PAH-emitting 
normal starbursts, and (2) energetically dominant, highly dust-obscured, cen- 
trally concentrated energy sources with no PAH emission. We interpret the latter 
component as a buried AGN. The fraction of ULIRGs showing some buried AGN 
signatures is higher in LINER ULIRGs than in Hll-region ULIRGs. Most of the 
luminous buried AGN candidates are found in ULIRGs with cool far-infrared 
colors. Where the absorption-corrected intrinsic AGN luminosities are derivable 
with little uncertainty, they are found to be of the order of lO^^L©, accounting for 
the bulk of the ULIRGs' luminosities. The 5-35 /im spectroscopic starburst/AGN 
classifications are generally consistent with our previous classifications based on 
3-4 nm spectra for the same sample. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: nuclei — galaxies: 
Seyfert — galaxies: starburst — infrared: galaxies 



Introduction 



Galaxies which radiate very large luminosities as infrared dust emission (Lir > IO^^Lq) 
are called ultraluminous infrared galaxies (ULIRGs; Sanders & Mirabel 1996). The huge 
infrared luminosities mean that: (1) very powerful energy sources, starbursts and/or active 
galactic nuclei (AGNs), are present hidden behind dust; (2) energetic radiation from the 
energy sources is absorbed by the surrounding dust; and (3) the heated dust grains re- 
emit this energy as infrared thermal radiation. In the local univer se (z < 0.3), the U LIRG 
population dominates the bright end of the luminosity function (jSoifer et al.l 119871 ). The 
importance of the ULIRG population i.e. the comoving infrared energy density, increases 
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rapidly with increasing redshift (ILe Floc'h et al.ll2005l ) . Identifying the dust obscured energy 
sources of nearby ULIRGs is thus very important, not onfy to unveil the nature of nearby 
ULIRGs, but also to clarify the obscured AGN-starburst connections in the high redshift {z > 
0.5) ULIRG population, where detailed investigations are difficult with existing facilities. 

In nearby ULIRGs, it is now widely accepted that the dominant energy sources are com- 
pact (<500 pc), highly dust-obscured nuclear c ores, rather than extended (>kp c), weakly- 
obscured starburst activity in the host galaxies (ISoifer et al.ll2000l : iFischen 120001 ). The most 
important issue is to determine whether the obscured compact cores are powered by very 
compact starbursts and/or AGNs, from which strong emission is produced in very compact 
regions around the central accreting supermassive blackholes. 

If a luminous AGN is present hidden behind dust in a torus geometry, as inferred for 
the classical Seyfert population (e.g., Antonucci 1993), a large amount of energetic AGN 
radiation can escape along the direction perpendicular to the torus, and form the so-called 
narrow line regions (NLRs; Antonucci 1993). NLRs are the sources of strong forbidden 
emission lines seen in the optical-infrared, and their relative strengths are different from 
those of a starburst because the spectrum at the ionizing high energy UV range (A < 912 
A) is harder in an AGN than in a starburst. An AGN obscured by a torus-shaped dust 
distribution, with a well-developed NLR, is thus distinguishable from a starburst galaxy 
relatively easily, through optical spectroscopy ( classified as Seyfert 2; Veilleux fc Oste rbrock 
1987) or high-resolution infrared spectroscopy (iGenzel et al.lll998l : lArmus et al.ll2007l ). 



Since the gas and dust in an AGN have angular momentum with respect to the central 
supermassibe blackhole, an axisymmetric spatial distribution is more natural than spherical 
geometry. In this case, the column densities can be high in certain directions (the torus 
directions), but low in others (along the torus axis). The classical Seyfert-type AGN is 
naturally understood as the population where most of the torus axis direction is transparent 
to the AGN's ionizing UV-soft-X-ray radiation. However, compared to the classical Seyfert 
AGNs, ULIRGs are knowii to co ntain a substantially larger amount of gas and dust in their 
nuclei ( ISanders fc Mirabel! Il996l ). For a fixed angular momentum, the increased amount of 
gas and dust results in higher column densities in all directions. If high density gas/dust 
blocks the bulk of the AGN's ionizing radiation at very small radii (<<10 pc) even in the 
torus axis direction, then no significant NLR will be present. It is very likely that any 
luminous AGNs in many ULIRGs exist in this buried condition, with no obvious Seyfert 
signatures in the optical spectra (hereafter "optical non-Seyfert"). Studying only AGNs 
with well-developed NLRs provides a very incomplete picture of the AGN population, and 
will miss many buried AGNs, whose understanding may be essential to clarify the true nature 
of ULIRGs. 
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An often-used indicator to separate a buried AGN from a normal starburst is the emis- 
sion from polycyclic a romatic hydrocarbons (PAHs), seen in infrare d spectra at Arcst = 3-25 
/xm in the rest-frame (iGenzel et aLlll998l : llmanishi fc Dudleyl |2000| ) . In a normal starburst, 
UV emitting Hll-regions and molecular gas and dust are spatially well mixed. PAHs are 
excited by far-UV photons from stars, and strong PAH emission is produced ir i photo- 
disso ciation regions (PDRs), the interface s between Hll-region s and molecular gas (jSellgren 
19811 ). unless the metallicity is very low (IHouck et al.ll2004bl ). In a normal starburst, the 
equivalent width of the PAH emission is insensitive to level of dust extinction, as both the 
PAH emission and the nearby continuum are similarly attenuated. Thus, a normal starburst 
should always show large equivalent width PAH emission, regardless of dust extinction 



Near an AGN, the PAHs are destroyed by strong X-ray radiation from the AGN (jVoitlll992 



Siebenmorgen et al.l 120041 ). If the PAHs are sufficiently shielded by obscuration from the X- 
ray emission of the central AGN, the PAHs can survive. However, the PAH-exciting far UV 
emission from the AGN will also be attenuated by the obscuration, so in a pure AGN with- 
out starburst activity, the PAH emission will be very weak. Instead, a PAH-free continuum 
from hot, submicron-sized dust grains heated by the AGN is observed. In a starburst/AGN 
composite galaxy, PAH emission from the starburst regions is seen, but its equivalent width 
will be smaller than in a pure starburst, because of the dilution by PAH-free continuum pro- 
duced by the AGN. Thus, this PAH-based method can be used to disentangle a buried AGN 
from a normal starburst at an obscured ULIRG's core. Since the PAH emission features are 
spectrally very broad, low-resolution infrared spectroscopy is adequate. In addition, ULIRGs 
are by definition bright in the infrared, so that this infrared PAH technique can be applied 
to a large number of ULIRGs in a reasonable amount of telescope time, and thus can be a 
very effective tool to systematically identify the hidden energy sources of ULIRGs' cores. 

Using low-resolution infrared spectroscopy, we also have another powerful tool to dis- 
tinguish buried AGNs from normal starbursts in the cores of ULIRGs, by investigating the 
geometry between energy sources and the gas and dust. In a n ormal starburst, the energy 



sources (stars) and gas/dust are spa tially well mixed (Figure la) (|Puxleylll99ll : iMcLeod et al. 



19931 : iForster Schreiber et al.ll200ll ). while in a buried AGN, the energy source (the central 
accreting supermassive blackhole) is very con apact and is more centrally concentrated than 



the surrounding gas and dust (Figure lb) (ISoifer et al.l l2000l : llmanishi fc Maloneyl 12003 



Siebenmorgen et al.ll2004l ). The difference in geometry is reflected in two different ways in 
the observed low-resolution infrared spectra. The first is that while the absolute optical 
depths of dust absorption features in the 3-10 fim range cannot exceed a certain threshold 
in a normal starburst with mixed dust /s ource geometry (Figure la) , they can be arbitrar- 
ily large in a buried AGN (Figure lb) (llmanishi fc Maloneyl l2003l : llmanishi et al.l l2006al : 
Levenson et al.l 120071 ). The second is that a strong dust temperature gradient is found in a 
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buried AGN (Figure 2), but not in a normal starburst. This dust temperature gradient can 
be investigat ed by comparing optical depths of dust absorption features at different infrared 
wavelengths ( Dudley fc Wynn- Williams! 119971 : llmanishil l2000l : llmanishi et al.l l2006al ) . 



In this paper we present low-resolution 5-35 am spectra of nearby ULIRGs at z < 0.15, 
to search for the signatures of optically elusive (IMaiolino et al.l 120031 ) buried AGNs in the 
cores of ULIRGs and to determine the energetic contribution of buried AGNs. Throughout 
this paper, Hq = 75 km s~^ Mpc~^, Qm = 0.3, and Q\ = 0.7 are adopted. The physical scale 
of 1" is 0.34 kpc in the nearest source at 2; = 0.018, 1.72 kpc at z = 0.1, and 2.44 kpc in the 
furtherest source at z = 0.15. 



2. Targets 



Our main targets are op tically non-Seyfert ULIRGs at 2; < 0.15 in the IRAS 1 Jy 
sample (IKim fc Sanders! 119981 ). In this 1 Jy sample, there are 69 ULIRGs at z < 0.15. 
Based on the optical classifications by Veilleux et al. (1999a; their Table 2), 28 ULIRGs 
are classified optically as LINERs and 20 ULIRGs are classified optically as HII regions. 
These 48 (=28+20) ULIRGs are selected, comprising ~70% of the 1 Jy sample at z < 0.15. 
The remaining ~30% U LIRGs are classified o ptically as Seyferts (17 sources) and optically 
unclassified (4 sources) (IVeilleux et al.! !l999al ). Table 1 summarizes the basic information 
and IRAS-hased infrared properties of the selected LINER and Hll-region ULIRGs. Low- 
resolution infrared 5-35 /im spectra for this complete sample are presented in this paper. For 



most s ources, ancillary ground-based infrared 3-4 /xm spectra are available (llmanishi et al. 
2006ah . 



3. Observations and Data Analysis 



O bservations of all th e 48 ULIRGs were performed usi ng the Infrared Spe ctrograph 
(IRS) (IHouck et al.!!2004al ) onboard Spitzer Space Telescope (jWerner et al.l!2004l ). All four 
modules, Short-Low 2 (SL2; 5.2-7.7 /im) and 1 (SLl; 7.4-14.5 /im), Long-Low 2 (LL2; 14.0- 
21.3 /xm) and 1 (LLl; 19.5-38.0 /zm) were used, to obtain full 5-35 nvd low-resolution (R ~ 
100) spectra. 13 out of the 48 sources were in the GTO target list (PI. J. Houck, pid = 105), 
and their archival spectra were analyzed. The remaining 35 sources were observed through 
our GO-1 program (PI = M. Imanishi, pid = 2306) and another GO-1 program (PI = S. 
Veilleux, pid = 3187). The GO-1 program by Dr. Veilleux covered only SL2 and SLl. For 
these ULIRGs, LL2 and LLl spectra were taken in our GO-1 program. 
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For our program, each observation consists of 60 sec x 2 cycles for SL and 30 sec x 
4 cycles for LL. The observation was performed both at A and B beam positions, and so 
total exposure time was 240 sec for both SL and LL. For the programs by Drs. Houck and 
Veilleux, the total exposure times are similar or smaller for bright sources. Table 2 tabulates 
the observing log of the 48 ULIRGs. The slit width is 3'.'6 or 2 pixels for SL2 (I'.'S pixel"^) 
and is 3'.'7 or ~2 pixels for SLl (I'.'S pixel"^). For LL2 and LLl, the slit widths are lO'.'S and 
10'.'7, respectively, which correspond to ~2 pixels for both LL2 (S'.'l pixel"^) and LLl (S'.'l 
pixel"^). 

The latest pipeline-processed data products at the time of our data analysis (Sll-14, 
pbcd files) were used. Frames taken at the A position were subtracted from those taken at 
the B position, to remove background emission, mostly the zodiacal light. Then, spectra were 
extracted in a standard manner using the task "apall" in IRAF 0. Apertures with 4-6 and 
4-5 pixels were employed for SL and LL data, respectively, depending on the spatial extent 
of individual sources. Then, spectra extracted for the A and B positions were summed. 
Wavelength calibration was made based on the files of the Spitzer pipeline processed data, 
named "bO_wavsamp.tbl" and "b2_wavsamp.tbl" for SL and LL, respectively. It is believed 
to be accurate within 0.1 /zm. A small level of possible error in wavelength calibration will 
not affect our main discussions. Since emission from all ULIRGs is dominated by spatially 
compact sources at the observed wavelength, flux calibration was performed using the Spitzer 
pipeline processed files "bOJluxcon.tbl" (SL) and "b2_fluxcon.tbl" (LL). For SLl spectra, 
data at Aobs > 14.5 fim in the observed frame are invalid (Infrared Spectrograph Data 
Handbook Version 1.0) and so were discarded. For LLl spectra, we used only data at Aobs 
< 35 /im, because the scatter is large at Aobs > 35 /im. No defringing was attempted. 

For flux-calibration, we adopted the values of the pipeline processed data. We made no 
attempt to re-calibrate our spectra using /i?^ 5" measurements at 12 /im and 25 /im, primarily 
because (1) only upper limits are provided for IRAS 12 /im photometry in a sizable fraction 
(>60%; 31/48) of the observed ULIRGs and for 25 fim data in some fraction (~20%; 10/48) 
of them (Table 1), and (2) the vahdity of re-calibrating IRS fixed-aperture slit spectra, using 
the large aperture photometry of IRAS, is not obvious for ULIRGs which possibly show 
a non-negligible fraction of spatially extended dim emission outside the IRS slit. Hence, 
the absolute flux calibration is dependent on the accuracy of the pipehne processed data, 
which is taken to be <20% for SL and LL (Infrared Spectrograph Data Handbook). This 
level of possible flux uncertainty will not affect our main discussions significantly. In fact, in 



^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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the majority of ULIRGs with IRAS detection at 25 fim, the Spitzer IRS 25 fim flux agrees 
within 30% to the IRAS 25 fim data, and the Spitzer IRS flux is generally smaller. In a 
small fraction of ULIRGs, Spitzer IRS 25 fim flux is ~50% of IRAS 25 /im data, possibly 
because extended dim emission is important. For ULIRGs with IRAS non-detection at 25 
/xm, the measured Spizter IRS flux at 25 /im is always smaller than the upper limit of the 
IRAS data. 

For a fraction of the ULIRGs, slight flux discrepancies between SLl and LL2 were 
discernible, ranging from 10% to 40%. When the discrepancy is seen, the SLl flux (3'.'7 wide 
slit) is usually smaller than the LL2 flux (10'.'5), suggesting the presence of extended emission 
in these ULIRGs. In these cases we adjusted the smaller SLl (and SL2) flux to match the 
larger LL2 flux. 

The main part of the 9.7 fim silicate absorption feature extends Arest = 8-13 /xm in the 
rest wavelength (e.g., Chiar & Tielens 2006). For ULIRGs at 2; > 0.11, the feature is not 
fully covered with SLl spectra at Aobs = 7.4-14.5 /im. Even for ULIRGs at z < 0.11, the 
continuum data at Arest = 13-14 /im, which are required to estimate the strength of the 9.7 
/im silicate absorption, fall beyond Aobs > 14.5 /im, outside the wavelength range of SLl. 
Thus, for the majority of the observed ULIRGs, LL2 spectra are essential to measure the 
optical depth of the 9.7 /im silicate absorption feature with a small ambiguity. 

Appropriate spectral binning with 2 or 4 pixels was applied to reduce the scatter of 
data points at SL2 (5.2-7.7 /im) for some faint ULIRGs, and at Aobs ~ 10 /im for ULIRGs 
displaying very strong 9.7 /im silicate absorption features. 



Result 



Figures 3 and 4 present low-resolution 5-35 /im spectra of ULIRGs classified optically 
as LINERs and Hll-regions, respectively. For most of the sources, full 5-35 /im spectra are 
shown here for the first time. Although spectra of the five ULIRGs in the GTO program 
(Arp 220, IRAS 08572+3915 12112+0305 1434 8-1447, and 22491-1808) were presented 
elsewhere (ISpoon et al.l l2006l : lArmus et al.l 120071 ) , we include their spectra to comprise a 
complete sample, to which our detailed and uniform analysis is applied. For these ULIRGs, 
their and our spectra are, in overall, consistent to each other. 

The spectra in Figures 3 and 4 are suitable for displaying the properties of the 9.7 
/im and 18 /im silicate dust absorption features. However, they are not very useful for PAH 
emission features. Figure 5 presents zoom-in spectra at Aobs = 5.2-14.5 /im of optical LINER 
and Hll-region ULIRGs, to better examine the properties of the PAH emission features. 
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4.1. PAH emission 



The majority of ULIRGs in Figure 5 show clearly detectable PAH emission features 
at Arest = 6.2 fim, 7.7 fim, and 11.3 /im. A fraction of the ULIRGs also display the PAH 
emission feature at Arest = 8.6 /xm. Emission at Arest = 12.8 /xm is also detected (Figures 3 
and 4), but it can be attributed to both the 12.8 /im PAH and [NeH] emission lines, which 
are not clearly resolved in our low-resolution spectra. We thus will not discuss the 12.8 /xm 
PAH emission feature further. 

Th e 6.2 Atm PAH emission feature is attributed to an aromatic C-C stretching funda- 
mental (lAUamandola et al.lll989l ). It is relatively weak, but is isolated both from other PAH 
emission and from the strong 9.7 fim silicate dust absorption feature. Although absorption 
at 6.0 /i m due to H9.O ice and at 6.85 /im by hydrogenated amorphous carbon (HAG) are 
present ( ISpoon et al.l 120021 ). their relatively weak nature, compared to the 9.7 /im silicate 
absorption, makes it possible to measure the 6.2 /im PAH emission strength with little ambi- 
guity. For this reason, the 6.2 /im feature has often bee n used to measure the PAH emission 
strength in a ULIRG (jFischerll2000l : ISpoon et al.ll2007l ). To estimate the strength of the 6.2 
/im PAH emission feature, we adopt a linear continuum determined from data at Arest = 6.1 
/im and 6.45 /im, which is shown as solid lines in Figure 5. 

The 11.3 /im PAH e mission is due to an out-of-plane aromatic G-H bending vibration 
(lAUamandola et al.lll989l ). Since it is situated close to the peak of the strong 9.7 /im silicate 
absorption, its flux attenuation can be severe if the PAH-emitting regions are highly dust 
obscured. Thus, its flux relative to other PAH features can be a good indicator of the dust 
obscuration of the PAH-emitting regions. To estimate the 11.3 /im PAH strength, we adopt 
a linear continuum connecting data points at Arest = H-O /im and 11.6 /im, which is shown 
as solid lines in Figure 5. 



The 7.7 /im PAH feature is due to aromatic G-G stretching vibrations (lAUamandola et al. 



19891 ). Although it is the strongest PAH emission feature, its substantial wavelength over- 
lap with the broad, strong 9.7 /im silicate dust absorption feature makes it very difficult 
to distinguish between 7.7 /im PAH emission and 9.7 /im silicate absorption, particularly 
for sources showing strong silicate absorption, such as the majority of ULIRGs. To reduce 
the effects of the 9.7 /im absorption feature, we adopt an excess above a linear continuum 
determined from data at Arest = 7.3 and 8.1 /im (shown as solid lines in Figure 5) as the 
7.7 /im PAH emission. Since many previous papers (e.g., Genzel et al. 1998) measured the 
7.7 fim PAH fluxes based on different definitions, readers must be careful when our value is 
compared with other estimates in the literature. In this paper, although the 7.7 /im PAH 
strength is shown for reference, it will not be used when possible uncertainty coming from 
the effects of the strong 9.7 /im silicate feature may complicate quantitative discussions. 
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The 8.6 /iin PAH emission due to in-plane aromatic C-H bending vibration (lAUamandola et al 



19891 ) is not used in this paper, because it is weak and overlaps with the 9.7 fim silicate ab- 
sorption, which make its strength highly uncertain. Other weak PAH emission features may 
be present (e.g., Arest = 5.7 /xm; Spoon et al. 2002), but are not discussed in detail because 
no new meaningful information is expected to be obtained. 

In summary, the 6.2 fim and 11.3 fim PAH emission features are the main focus of 
our discussion, and the strength of the 7.7 /im PAH feature based on our own definition is 
shown just for reference. We fit the PAH emission features with Gaussian profiles, which 
reproduce the observed data reasonably well. The observed rest-frame equivalent widths and 
luminosities of the 6.2 /xm, 7.7 /xm, and 11.3 /im PAH emission features, based on our adopted 
continuum levels, are summarized in Table 3. The uncertainties coming from the fittings are 
unlikely to exceed 30%, so once the continuum levels are fixed, the total uncertainties are of 
this order. 



4.2. Silicate absorption 

Based on our adopted continuum levels, the 9.7 /im silicate absorption feature is found 
in almost all ULIRGs, and 18 fim silicate absorption is also seen, particularly in ULIRGs with 
strong 9.7 /xm absorption (Figures 3 and 4). Silicate dust grains show two strong absorption 
peaks centered at Arest = 9.7 fim and 18. 5 fim. Outside these p eaks, a weak absorption tail 
extends from Arest = 8 /im to >30 /im Jchiar fc TielensI l2006h (see also Figure 6a in this 



paper). Thus, all wavelengths between Arest = 8-30 /im show excess absorption, compared 
to the extrapolation of dust extinction curve from wavelengths outside Arest = 8-30 /im. 
Hence, to accurately estimate the optical depths of the 9.7 /im (rg.y) and 18 /xm (rig) silicate 
absorption features, a continuum level must be determined using data at Arest < 8 /im and 
>30 /im. At Arest = 6-8 /im, PAH emission can contaminate the continuum level, so that 
in practice data points at Arest ~ 5.6 /im are used at the shorter wavelength side. Data 
points at Arest. ^ 7.1 /im are also less affected by PAH emission and absorption features 



(jSpoon et al.ll2006l ). and so can be used as additional continuum points. Table 4 summarizes 
the optical depths of the 9.7 /xm silicate absorption (rg.y), against power-law continuum levels 
determined from data points at Arest = 5.6 /im, 7.1 /im, and Aobs = 34-35 /im (the longest part 
of the IRS spectra) shown as the dashed lines in Figures 3 and 4. A power-law continuum is 
assumed here, because the intrinsic unabsorbed 5-35 (j, m spectrum of AGN-he ated hot dust 



emission is roughly approximated by a power-law form (jWeedman et al.ll2005l ). 



The wavelengths of these data points used for the continuum determination are, how- 
ever, far away from the 9.7 /im and 18 /im silicate peaks. The optical depth values of the 
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9.7 /iin and 18 silicate absorption features can vary depending on the assumptions of 
the continuum data points as well as the continuum shape (power-law, linear, spline, or 
polynomial etc). Additionally, the power-low continua do not necessarily provide perfect fits 
to observed data points outside the silicate features. In the spectrum of the weakly obscured 
ULIRG IRAS 01004-2237 (Figure 4, third plot), data points at Aobs = 13-18 /im are clearly 
above the power-law continuum (dashed line), indicating that the power-law continuum only 
roughly represents the overall continuum. The relative ambiguity of the rg ^ value is small 
because the 9.7 fim silicate feature is intrinsically strong. However, the optical depth of 
the much weaker 18 /im silicate absorption feature can be highly uncertain. Consequently, 
discussions of the presence or absence of strong dust temperature gradients, based on the 
relative strength of the 9.7 fim and 18 fim silicate dust absorption features (Figure 2), could 
also be highly uncertain, if we use rg 7 and rig as defined above. 

For this reason, we define Tq j and r{g. The Tg 7 value is the optical depth of the 9.7 fim 
silicate feature against a power-law continuum determined from data points at Arest = 7.1 
fim and 14.2 /im. The r{g value is the optical depth of the 18 /xm silicate feature against 
a power-law continuum determined from data points at Arest = 14.2 /im and 24 /im. These 
continua are shown as dotted lines in Figures 3 and 4. With these new definitions, the 
ambiguities are reduced, because the continuum levels are determined using data points just 
outside the 9.7 /im and 18 /im features, closer to the absorption peaks. The Tg 7 values are 
shown for all ULIRGs in Table 4 (column 3). The r{g values are also shown in Table 4 
(column 4) for ULIRGs where the 18 /im feature is clearly seen in absorption. 

The r^g values are appropriate to estimate the strength of the 18 /im silica te feature, if 



the 1 8 /zm profile is similar to that observed in the ULIRG IRAS 08572-1-3915 (jSpoon et al. 



20061 ) . However, the 18 /xm silicate profile of the Galactic, highly obscured star GCS3- I 



extends far beyond Arest ~ 24 /im, and reaches to Arest ~ 29 /im (jChiar fc Tielend l2006l ) 
(see also Figure 6a in this paper). Thus, we also define r(g, which is the optical depth of 
the 18 /im feature against a power-law continuum determined from data points at Arest = 
14.2 /im and 29 /im (shown as dashed-dotted lines in Figure 3 and 4). The r"g values are 
also summarized in Table 4 (column 5) for ULIRGs with clearly detectable 18 /im silicate 
absorption. 



4.3. Ice absorption 

A significant fraction of ULIRGs display dips at the shorter wavelength part of the 6.2 
/im PAH emission features. To better visualize the dips. Figure 7 presents zoom-in spectra 
at Aobs = 5.2-9 /im for ULIRGs showing clear or possible dips. Spectra of two ULIRGs 
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without discernible dips (IRAS 17044+6720 and 00456—2904) are also shown as examples 
of no clear detection. We ascribe the dips to the 6.0 fim H^ O absorption feature (bending 



mode ), as has been found in both Galactic obs cured stars (iGibb et al.ll2000l : iKeane et al. 



200ll ) and obscured ULIRGs (ISpoon et al.ll2002l ). Since the 6.0 /im ice absorption feature 
is sandwiched by PAH emission features at Arest = 5.7 and 6.2 /im, it is sometimes difficult 
to distinguish whether the dip is due to ice absorption or to PAH emission features. For 
ULIRGs with clearly detectable 6.0 /im H2O ice absorption features, their observed optical 
depths (Tg.o) are summarized in Table 5. 

For ULIRGs showing clear ice (and HAG) absorption features, the continuum flux levels 
just outside of the 6.2 /xm PAH emission are reduced. For these ULIRGs, approximate 
EW6.2PAH values after these corrections are also shown in Table 3. 



Discussion 



5.1. Magnitudes of detected starbursts 

PAH emission luminosities can roughly probe the absolute magnitudes of PAH-emitting 
normal starbursts, if the extinction of the starbursts is say less than about 20-30 mag in 
Av, because this amount of Ay produces little flux atte nuation (<1 mag) at Ar -e st > 5 Atm, 
aside from the strong 9.7 fim silicate absorption feature (IRieke &: Lebofskylll985l : iLutz et al. 
1996h . 



Table 3 (columns 9 and 10) tabulates the 6.2 fim PAH to infrared luminosity ratio, 
L6.2PAh/Lir, and the 11.3 /im PAH to infrared luminosity ratio, Lh.spah/Lir. The ratios in 
normal starburst galaxies wi th modest dust obsc uration (Ay < 20-30 mag) are estimated to 
be L^.2pah/Lir ~ 3.4 x 10"^ (jPeeters et al.i2004l ) and Lii.3pah/Lir ~ 1.4 x 10"^ (jSoifer et al. 



2OO2I ) . Figure 8a and 8b display the PAH-to-infrared luminosity ratio (ordinate) and PAH 
equivalent width (abscissa) for the 6.2 fim and 11.3 /im PAH features, respectively. The 
observed L6.2pah/Lir and Lh.spah/Lir ratios in the majority of the observed ULIRGs are 
systematically lower than those seen in modestly-obscured normal starburst galaxies, sug- 
gesting that the detected starburst activity as measured by the PAH emission in ULIRGs is 
energetically insignificant. The dominant energy sources of ULIRGs must therefore be either 
(1) highly-obscured (Ay >> 30 mag) starbursts, where the emitted PAH luminosities are 
substantially attenuated by dust extinction, or (2) buried AGNs, which produce intrinsically 
weak PAH emission. 
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5.2. Candidates of ULIRGs possessing powerful buried AGNs 

5.2.1. Low equivalent widths of PAH emission features 

Whether the dominant energy sources of ULIRGs' cores are highly-obscured normal 
starbursts with mixed dust/source geometry (Figure la), or buried AGN with centrally- 
concentrated energy source geometry (Figure lb,c), can be distinguished using the equivalent 
width of the PAH emission. Since the PAH equivalent width (EWpah) must always be large 
in a normal starburst regardless of the amount of dust extinction, a small EWpah value 
requires a contribution from a PAH-free continuum-emitting energy source (§1). 

Figure 9 (a), (b), and (c) compare the equivalent widths of the 6.2 /im, 7.7 /im, and 
11.3 /im PAH emission features. Positive correlations are found in all plots, suggesting that 

(1) variation of the relative PAH emission strength is modest, with some level of scatter, and 

(2) PAH-emitting starbursts usually show strong PAH features in all the emission modes at 
6.2 /xm, 7.7 /xm, and 11.3 /xm. 



Brandl et al.l (120061 ) have investigated the 6.2 /xm and 11.3 /xm PAH emission equiva- 
lent widths, based on similar definitions of PAH strength to ours, in nearby, normal star- 
burst galaxies without obvious Seyfert signatures in the optical (18 sources), and found that 
EW6.2PAH = 561±73 nm (median 553 nm) and EWh.spah = 617±187 nm (median 597 nm). 
The scatter is larger for the 11.3 /im PAH emission feature than the 6.2 /im PAH feature. 
No starburst galaxies show values smaller than EW6,2pah = 459 nm and EWh.spah = 316 
nm. Since the PAH equivalent widths are indep endent of the infra red luminosities of the 



starburst galaxies in the range Ljr =10 (iBrandl et al.ll2006l ). we apply these values 

to our ULIRG sample (Lir > lO^^Lg 



The median values for the 6.2 /im PAH emission equivalent widths (rest-frame) are 
EWfi 2PAH = 230 nm for the whole of our ULIRG sample (LINER + Hll-region), 225 nm for 
LINER ULIRGs, and 300 nm for Hll-region ULIRGs. The median 11.3 /xm PAH equivalent 
widths are EWh.spah = 390 nm (LINER + Hll-region), 350 nm (LINER) and 460 nm (HH- 
region). The median equivalent widths o f both PAH features are smaller in ULIRGs than 



in starburst galaxies (IBrandl et al.l 120061 ). suggesting that PAH-free continuum emission, 
originating in AGN- heated hot dust, contributes more strongly to the observed fiuxes in 
ULIRGs. We assume that if a PAH equivalent width is more than a factor of ~3 smaller than 
the median value in normal starburst galaxies, then an AGN-originating PAH-free continuum 
contributes importantly to the observed 5-12 /xm fiux and dilutes the PAH emission. With 
this criterion, ULIRGs with EW6.2pah < 180 nm and EWh.spah < 200 nm are classified as 
sources displaying clear signatures of luminous AGNs. Since these adopted values are clearly 
below the minimum values observed in normal starburst galaxies, we believe the criterion is 
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reasonable and conservative to pick up luminous buried AGN candidates. 

For the 7.7 iim PAH equivalent widths, based on our definition of the PAH strength 
(§4.1), the median 7.7 iim PAH equivalent widths in our ULIRG sample are EWy.ypAH = 490 
nm (LINER) and 695 nm (HIT- region). In our ULIRG sample, the median PAH equivalent 
widths are smaller in LINER ULIRGs than in Hll-region ULIRGs for all of the 6.2 //m, 
7.7 /xm, and 11.3 //m features. We classify ULIRGs with EW7.7PAH < 230 nm, one-third of 
the median value for Hll-region ULIRGs, as those showing luminous buried AGN signatures. 
Since some fraction of Hll-region ULIRGs may contain luminous buried AGNs, this EW7.7PAH 
threshold is even more conservative than those for the 6.2 /xm and 11.3 /xm PAH emission 
features, as determined above. 

Table 6 (columns 2-4) presents detection or non-detection of buried AGN signatures 
based on the PAH equivalent width threshold. The fraction of ULIRGs with clearly de- 
tectable buried AGN signatures is 16/48 (33%; LINER + Hll-region), 10/28 (36%; LINER) 
and 6/20 (30%; HH-region) based on observed EW6.2PAH, or 20/48 (42%; LINER + HII- 
region), 13/28 (46%; LINER) and 7/20 (35%; HH-region) based on EW6.2pah after correction 
for the ice and HAG absorption feature. The fraction is 4/48 (8%; LINER + HH-region), 
3/28 (11%; LINER) and 1/20 (5%; HH-region) based on observed EW7.7PAH, 12/48 (25%; 
LINER + Hll-region), 8/28 (29%; LINER) and 4/20 (20%; HH-region) based on observed 
EWii.3PAH- The fraction of ULIRGs where very low PAH equivalent widths are found in 
any of the features is 19/48 (40%; LINER + HH-region), 13/28 (46%; LINER) and 6/20 
(30%; HH-region). If the corrected EW6.2PAH values arc adopted, it is 22/48 (46%; LINER 
+ HH-region), 15/28 (54%; LINER) and 7/20 (35%; HH-region). 

Since we adopted a conservative threshold, this method provides only candidate ULIRGs 
where PAH-free continuum emission from buried AGNs contribute very importantly to the 
observed IRS spectra. In starburst/buried- AGN composite ULIRGs (Figure Ic), PAH emis- 
sion from starbursts surrounding the central buried AGNs is less obscured and less flux- 
attenuated, making the observed PAH equivalent widths relatively large, particularly when 
AGNs are deeply buried. Such buried AGNs should be searched for using other methods. 

5.2.2. Absolute optical depths of dust absorption features 

In a normal starburst with mixed dust/source geometry (Figure la), the foreground, 
less-obscured, less-attenuated emission (which shows only weak dust absorption features) 
dominates the observed flux, and therefore the observed optical depths of dust absorption 
features cannot exceed a certain threshold, unless very unusual dust composition patterns are 
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assumed ( Imanishi fc Maloney 20031 : Imanishi et al. I l2006ah . In a buried AGN, with centrally- 
concentrated source geometry (Figure lb), the so-called foreground screen dust model is 
applicable, and the observed optical depths can be arbitrarily large. Hence, detection of 
strong dust absorption features, with optical depths substantially larger than the upper 
limit set by the mixed dust /sour ce geometry, argues for a foreground screen dust geometry, 
as expected from a buried AGN (jimanishi fc Maloneyll2003l : llmanishi et al.ll2006al ). 



To differentiate the geometry using the optical depths of dust absorption features, the 
9.7 fim silicate dust absorption is particularly useful, because (1) it is much stronger than the 
18 fim sihcate feature (Figure 6a), and (2) its in t rinsic streng th has been widely investigated 
in Galactic obscured objects (IRoche &: AitkenI Il984l . llQSSi ). Assuming the Galactic dust 
extinction curve, the upper limit for the optical depth of the 9.7 /im silicate absorption 
feature in the mixed dust/source geometry can be derived. Appendix A presents detailed 
calculations and we obtain a maximum value of Tg ^ < 1.7 for mixed dust /source geometry. 
In fact, none of the starburst ga laxies (Ltr < lO^^f Lp-,) lacking optical Seyfert signatures 
so far observed with Spitzer IRS (IBrandl et al.ll2006l ) show values that exceed the threshold 
(tq y = 1.7). Thus, it is reasonable to assume that ULIRGs with Tg j substantially larger 
than 1.7 possess buried AGNs with centrally-concentrated energy source geometry. 

In ULIRGs, it is possible that the dust extinction curve is very different from the Galac- 
tic one, because dust coagulation in the high density environment in ULIRGs' cores could 



increase typical dust g rain size up to a few microns (iLaor &: Draind Il993l : iMaiolino et al. 



2000al lbl: I Imanishi! I2OOII ). However, we believe this possibility has minimal effect on our con- 
clusions about the geometry, as long as we adopt Xg y (peak optical depth of the 9.7 /im 
feature, relative to the nearby continuum just outside the feature; see §4.2), rather than 
Tg 7 (peak optical depth of the 9.7 /im feature and continuum, relative to continuum outside 
8-30 /xm; see §4.2). Figure 6b shows the absorption optical depth of astronomical silicate 
dust grains of three different sizes. The absolute absorption optical depth at the peak of 
the 9.7 fim feature becomes stronger with increasing dust grain size. However, both the 
absorption optical depth inside the 9.7 fim feature and the nearby continuum just outside 
the feature increase in a very similar way with increasing dust size, so that the strength of 
9.7 /xm feature relative to the nearby continuum (rg^), is virtually unchanged with varying 
dust grain size. The argument of the geometry based on Tg y values are therefore robust to 
possible increase of dust grain size in ULIRGs' cores. We conservatively take ULIRGs with 
Tg 7 > 2 as candidates for harboring luminous centrally-concentrated buried AGNs. 

Individual ULIRGs with clearly detectable buried AGN signatures based on large Tg j 
values (>2; Table 4) are marked with open circles in Table 6 (column 5). The median Tg 7 
values are 1.8 (LINER + Hll-region), 2.1 (LINER), and 1.6 (Hll-region). The fraction of 
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ULIRGs with Tg^ > 2, indicative of buried AGNs with centrally-concentrated energy source 
geometry, is 20/48 (42%; LINER + Hll-region), 14/28 (50%; LINER), and 6/20 (30%; HII- 
region). The fraction of ULIRGs showing buried AGN signatures is again higher in LINER 
ULIRGs than Hll-region ULIRGs, as was previously found based on the PAH-equivalent- 
width method (§5.2.1). 

We must raise three cautions regarding this method. First, a Xg ^ value smaller than 1.7 
does not necessarily preclude the presence of a luminous AGN, because a weakly obscured 
AGN, of course, shows a small Tg 7 value. The best example is IRAS 01004—2237 (rg.7 ~ 0.4), 
which shows very low PAH emission equivalent widths, as expected from an AGN-dominated 
source. 

Next, we used the threshold of Tg ^ > 2, based on the observed values. For ULIRGs 
with detectable PAH emission from modestly obscured starbursts, the observed fluxes are 
a superposition of the emission from the highly obscured cores and the starburst emission 
(Figure Ic), and so the observed Tg j values are reduced compared to the Tg 7 values only 
toward the highly obscured energy sources at ULIRGs' cores. Since subtraction of the mod- 
estly obscured starburst component, using a starburst template spectrum, introduces an 
additional uncertainty, we do not attempt it. Thus, we may miss ULIRGs which show large 
Tg 7 values (>2) toward the cores but small observed Tg 7 values (<2) due to the starburst 
dilution. 

Finally, a starburst nucleus obscured by foreground absorbing dust (Figure Id) can also 
show a Tg 7 value larger than the threshold determined by the mixed dust/source geometry, 
because in this foreground screen dust model is applicable. Assuming Gala ctic 

dust properties, r9 7 > 2 corresponds to Ay > 20-4 mag (IRoche fc Aitkenlll984j . Il985l ). If 



ULIRGs are the mergers of gas-rich spiral galaxies (ISanders et al.lll988ar). gas and dust in 



the original galaxies can be qu ickly transferred to the nuclear regions (IBarnes fc Hernquist 



19961 : iMihos fc Hernquistlll996l ). The bulk of possibly newborn stars and resulting dust are 
also expected to distribute in the nuclear regions. In fact, most of high density gas, dust, 
and energy sources (starbursts and/or AGNs) in ULIRGs are observed to be concentrated to 



nuclea r cores with <kpc scale (IDownes &: Solomonlll998l : ISoifer et al.l I2OOOI : llmanishi et al. 



2006bl ). The amount of gas/dust in ULIRGs' hosts at >kpc scale (which are morphologically 
disturbed) is unlikely to substantially exceed that of the original merging gas/dust-rich 
galaxies. In our Ga laxy, a gas-rich sp iral, Ay value is estimated to be >30 mag in an 
equatorial direction (IRieke et al.lll989l ). and should be much smaller in other directions. 
It is thus not easy to produce such a large Ay value (>20-40 mag) with dust in a host 
galaxy, unless the host is viewed from a direction with particularly large dust column density. 
Although a small fraction of ULIRGs with rg7 > 2 may have this geometry (Figure Id), it 
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is very unlikely that the majority do. 



5.2.3. Strong dust temperature gradients 

A buried AGN with centrally-concentrated energy source geometry (Figure lb) should 
show a strong dust temperature gradient (Figure 2), whereas a normal starburst nucleus with 
mixed dust /source geometry does not, whether it is obscured by foreground dust (Figure 
Id) or not (Figure la). The presence of a strong dust temperature gradient thus sup- 
ports the buried AGN scenario, and it can be detected by comparing the optical depths of 



dust absorption features at different infrared wavelengths (IDudlev fc Wynn- Williams! Il997 



Imanishil 12000 : llmanishi fc Maloneyl l2003l : llmanishi et al.l l2006al ). As Figure 2 illustrates. 



in a buried AGN, the dust column density estimated from observations at Arest > 3 /zm 
probes deeper inside of the obscuring material as the wavelength is decreased {i.e., Av(3/im) 
> Av(10/im) > Av(20/im)}. In this geometry, dust with a fixed temperature range from a 
fixed layer contributes strongly to the observed continuum emission at each wavelength, with 
a small contribution from other components (Figure 2). Thus, in a starburst/buried AGN 
composite (Figure Ic), PAH-free continuum emission from the AGN- heated dust should be 
detectable even if its optical depth is larger than unity at a particular observed wavelength, 
provided the buried AGN is intrinsically sufficiently luminous that the flux-attenuated (e""^; 
r > 1) AGN-heated dust emission contributes signiflcantly to the observed flux. This has 
been demonstrated in the 3-4 /zm spectra of the ULIRGs IRAS 08572+3915 and UGC 5101 
{n-ium > 5; Imanishi et al. 2000; 2001; Imanishi & Maloney 2003). 

On the other hand, in a normal starburst (Figure la), no signiflcant dust temperature 
gradient is expected, and so the estimated dust extinction should be similar at different 
wavelengths {Av(3/im) ~ Av(10/im) ~ Av(20/im)}. Since in a normal starburst, emission 
arises uniformly within an emitting volume, we can typically probe emission with optical 
depth less than about unity. If the extinction is large enough that the short wavelengths 
do not probe the whole emitting region, the estimated dust extinction would decrease with 
decreasing wavelength {Av(3/im) < Av(10;um) < Av(20/im)}. Hence, the relations among 
Av(3/im), Av(10/xm), and Av(20/im) are totally different between a buried AGN (Figure lb) 
and a normal starburst (Figure la). 

For a buried AGN, the t[^/t'qj or t'I^/t^j ratio in the Spitzer IRS spectrum becomes 
small compared to the value in the absence of a strong dust temperature gradient. The 
'^is/'^Lr {'''is/'^d.i) ratio is also insensitive to possible changes in dust grain size in ULIRGs, 
because the relative strength of r{g (r"g) and Tq y varies little with varying dust grain size 
(Figure 6b). Furthermore, since both absorption features are due to silicate dust grains, the 
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ratio is unaffected by possibly different dust composition patterns in individual galaxies. 

For ULlRGs witli large contributions from modestly obscured PAH-emitting starbursts 
to observed spectra, the 18 /im silicate dust absorption features are diluted by the PAH 
complex at Arcst = 16.4 /im and 17.12 /x m, as well as H2 emission at Arcst = 17.04 /xm 



and [SIII] emission at Arcst = 18.71 /im (ISmith et al.l |2004| ). Subtraction of a template 



spectrum of the PAH-emitting starburst introduces an additional ambiguity in the estimate 
of the intrinsic r(g (rfg) value toward the ULIRG's core. To minimize the uncertainties, 
we investigate the r^g/rg y (Tfg/rg y) ratios only for ULIRGs where the contributions from 
PAH-emitting normal starbursts are small and the 18 fim silicate features are clearly seen in 
absorption in the observed spectra. The main aim of applying this T^g/rg 7 (T"g/rg 7) method 
to these ULIRGs is to further strengthen the case for buried AGN previously suggested by 
the large Tg j values, rather than to discover new buried AGN candidates. 

To investigate the presence of a strong dust temperature gradient in a ULIRG, based 
on a small observed t[q/tq y {tis/'^qj) ratio, we first have to know its intrinsic ratio in the 
absence of a temperature gradient. Figure 6a presents the observed profiles of two Galactic 



obscured stars, GCS3-I and WR98a (IChiar fc Tielensll2006l ). They are thought to be mostly 



obscured by foreground interstellar medium far away from the central illuminating stars, 
rather than circumstellar material, so no significant temperature gradient is expected. To 
conservatively argue that a strong dust temperature gradient is found, based on a small 
observed T^g/xg ^ (T^g/Tg y) ratio, we have adopted the profile showing the smallest intrinsic 
'^is/^Lt (''"is/'^g.?) ratio as a template. The intrinsic ratio of GCS3-1 is estimated to be t[^It!^ i 
= 0.3 and r^g/rg ^ = 0.35. The ratios of WR98a, and other Galactic obscured stars for which 
sihcate profiles are investigated in detail (OH-IR127.8+0.0, WR118, and WR112; Chiar & 
Tielens 2006), are larger than that of GCS3-I, so that adopting the observed GCS3-I profile 
as a template is appropriate here. If an observed T^g/rg^ or t^^/tqj ratio is substantially 
smaller than the adopted intrinsic value {t[q/tq j = 0.3 and t^s/tq^j = 0.35), then the most 
reasonable explanation is the presence of a strong dust temperature gradient. 

Table 4 (columns 6 and 7) summarizes the t[^/tq 7 and tI^/t!^ 7 ratios for ULIRGs show- 
ing clear 18 /im silicate absorption (mostly Tg 7 > 2). For a significant fraction of such 
ULIRGs, the observed ratios are smaller than the conservatively assumed intrinsic one, sug- 
gesting the presence of strong dust temperature gradients. This provides additional evidence 
for these ULIRGs for centrally-concentrated buried AGNs, previously suggested by low PAH 
equivalent widths and/or large r9 7 values (>2). 

ULIRGs showing detectable signatures of strong dust temperature gradients are marked 
as open circles in Table 6 (column 6). The fraction of ULIRGs with detectable signatures of 
strong dust temperature gradients is 7/48 (15%; LINER + Hll-region), 6/28 (21%; LINER), 



- 18 - 



and 1/20 (5%; HII- region). The fraction of ULIRGs with detectable buried AGN signatures 
is again higher in LINER ULIRGs than Hll-region ULIRGs. 

Again, there are several caveats. First, we adopted the observed profile of GCS3-I as 
the template, since the Tig/rgy ratio in GCS3-I is the lowest among five Galactic sources 
so far investigated in detail. A source with an even smaller intrinsic T^g/rg^ ratio may be 
discovered among Galactic sources in future observations. If ULIRGs showed intrinsically 
smaller T^g/xg j ratios than the adopted ratio from GCS3-I, then the small observed ratios in 
ULIRGs would not necessarily require the strong dust temperature gra dients. However, for 



IRAS 08572+3915, a bright ULIRG with a small r^g/r^ 7 ratio (Table 4). ISpoon et all fl2006f ) 
independently found the signature of a strong dust temperature gradient, based on the fact 
that the observed optical depths of crystalline silicate dust absorption at Arest ~ 11 /xm, 16 
/xm, 18.5 /im, 23 /im, and 28 /im are systematically weaker at the longer wavelengths than 
expected (their Figure 2 bottom). A similar trend was found also in the average spectrum 
of ULIRGs with Tg ^ > 2 (Spoon et al. 2006; their Figure 5). We thus believe that the 
small T^g/Tg 7 ratios observed in several ULIRGs are most naturally explained by strong dust 
temperature gradients. 

Next, although we investigated the T^g/rg^ ratios only for ULIRGs with small PAH 
contributions from modestly obscured starbursts, these small contributions could reduce both 
r(g and Xg 7 (§5.2.2), and slightly alter the t[^/t^^-j ratios, compared to a pure spectrum of 
the highly obscured energy sources in the cores. However, the average s pectrum of starburs t 



galaxies shows the 30 /^m to 5 /xm (rest-frame) flux ratio to be ~900 (iBrandl et al.ll2006l ). 
while the majority of the ULIRGs with small t[^/t!^^ show flux ratios larger than >900 
(Figures 1 and 2). That is, the average starburst spectral energy distribution is bluer at Arest 
= 5-30 /im than those ULIRGs with small t[^/t^'j- Assuming that the average spectrum 
is representative of the spectral energy distribution of the modestly obscured starbursts 
in these ULIRGs, this possible effect increases (not decreases) the observed t[^/t'q^j ratios, 
compared to an uncontaminated spectrum of the highly obscured core. Thus, the observed 
small T(g/r9 7 ratios in these ULIRGs are still best explained by strong dust temperature 
gradients around buried AGNs, with centrally concentrated energy source geometry. 

Finally, of course, silicate dust grains in the cores of ULIRGs may have intrinsically 
larger T^g/xg 7 ratios than the conservatively assumed one. In this case, if a ULIRG shows a 
decreased t[^/t!^,^ ratio compared to the intrinsic one, due to a dust temperature gradient, 
but still has an observed ratio similar to the assumed one, then such a ULIRG would not be 
classified as a buried AGN. 
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5.2.4- Comparison of the strength of silicate absorption and PAH emission features 

Figure 10a compares the optical depth of the 9.7 fim sihcate dust absorption (rg 7) and 
the rest-frame equivalent widths of the 6.2 fim PAH emission feature (EW6.2PAh)- With the 
exception of the two outliers with both small EW6.2pah and small Tg j (IRAS 01004—2237 and 
17044+6720, both weakly obscured AGNs), there is a weak anti-correlation, and ULIRGs 
with Tg^ > 2 (i.e., buried AGN candidates) tend to show small EW6.2pah compared to 
starburst-like ULIRGs with small Tg^ (<2). 

A highly-obscured (Ay >> 30 mag) normal starburst galaxy with mixed dust /source 
geometry (scenario A) predicts a relatively flat distribution (i.e., relatively constant EW6.2PAH 
values with varying Tg 7), because the PAH equivalent width is unchanged by dust extinction 
(§1). On the other hand, scenario B, where (1) the PAH emission comes from modestly 
obscured (Ay < 20-30 mag) starbursts exterior to the deeply obscured ULIRG cores, and (2) 
the large Tg 7 values reflect dust column density toward the buried AGNs in the cores, predicts 
an anti-correlation between Tg 7 and EW6.2pah- As the contribution from the AGN-produced 
PAH-free continuum to the observed flux increases, the the EW6.2PAH value decreases. At 
the same time, the increasing contribution from the buried AGN emission (with centrally- 
concentrated energy source geometry) results in the increase of the observed Tg 7 value, 
compared to a pure starburst with mixed source/dust geometry. The trend in Figure 10a 
(see also Spoon et al. 2007) supports the scenario B for ULIRGs with Tg 7 > 2. 

Figure 10b compares rg7 and the rest-frame equivalent widths of the 11.3 fim PAH 
emission feature (EWh.spah)- Compared to Figure 10a, the decreasing trend of EWh.spah 
at Tg 7 > 2 appears slight, at best. In scenario B, the flux depression of the AGN-produced 
continuum flux at Arcst = 11-3 /im is particularly large due to the strong 9.7 fim silicate 
absorption. However, the 11.3 /xm PAH emission flux is less affected by this absorption 
because the PAH emission comes from starburst regions outside the core. Consequently, the 
decrease of the EWh.spah value with increasing Xg 7 value is small compared to EW6.2pah, 
explaining the trends in both Figures 10a and 10b. 

Figure 11 compares the Tg 7 values with the 11.3 /im to 6.2 fim PAH luminosity ratio. In 
scenario A, the 11.3 /im to 6.2 /im f lux ratio should dec rease with increasing Xg 7, as actually 



seen in normal starburst galaxies (IBrandl et al.ll2006l ). because the flux depression of the 



11.3 /im PAH emission inside the strong 9.7 /im absorption is larger. On the other hand, in 
scenario B, no correlation is expected because both the 6.2 /im and 11.3 /im PAH emission 
come from the modestly obscured starbursts, unrelated to the observed Tg 7 values. Figure 
11 shows no significant correlation, supporting scenario B. 
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5.2.5. Is the buried AGN scenario unique': 



ULIRGs with large Tg 7 (>2) and small T^g/rg y ratios are naturally explained by the 
presence of buried AGN (scenario B) . We here briefly comment on two alternative scenarios 
which are possible but we believe are unlikely. 

First, in a normal starburst with mixed dust /source geometry, if each star is deeply 
buried in a large amount of circumstellar dust cocoon, such as a protostar, then a strong 
dust absorption feature (rg 7 > 2) could be produced. However, the time scale of such a 
protostar phase is only <10^ yr (jLadalll987l ). which is much less than the typical time scale 
of a normal starburst (~ 10''^^ yr; Moorwood 1996) or even the maximum starburst (~ lO'^^^ 
yr; Elmegreen 1999), or a ULIRG time scale (~ 10*^^^ yr; Murphy et al. 1996). It is therefore 
very unlikely that all dust-heating stars are protostars in a significant fraction of ULIRGs. 
Furthermore, this scenario produces no strong dust temperature gradient, and so cannot 
reduce the observed T^g/rg y ratio. 

Second, it may be argued that exceptionally centrally-concentrated starbursts remain 
a possibility (Figure le). However, while starburst galaxies with Lir =1O^ '^^^^'^L0 are well 
explained by the mixed dust/source geometry (Brandl et al. 2006; see §5.2.4 of our paper), 
there is no reason why starbursts at the cores of ULIRGs (Ljr > IO^^Lq) suddenly change 
their properties to exceptionally centrally-concentrated nature (Figure le). Furthermore, 
Soifer et al.l (I2OOOI ) studied the cores of seven nearby ULIRGs and found that the observed 
emission surface brightnesses in the majority (six out of seven) of these ULIRGs' cores are 
very high (~10^^ Lq kpc^^). Even if the energy sources are uniformly distributed over 
the cores (<500 pc in size), the observed surface brightnesses are c lose to the maximum 
values seen in the cores of Galactic HII regions (jSoifer et al. 2000). If an exceptionally 
centrally- concentrated starburst were the luminosity source in the core of a ULIRG, the 
emission surface brightness of the starburst would have to be extremely high, >10^^ L© 
kpc~^. Among known star-forming regions, such a high surface brightness could only barely 
be achieved with a sup er star cluster, whe re many young stars are formed within a area 
of less than a few pc (iGorjian et al.l I2OOII ). The absolute luminosity of such super star 
clusters is, however, only of order 1O^~^°L0, and so the super star cluster scenario may work 
for the centrally-concentrate d energy sources of low luminosity (Ljr ~ lO^'^Lo) galaxy nuclei 
( llmanishi &: Nakanishill2006l ). For the centrally-concentrated energy sources at ULIRG cores 
both high emission surface brightness and large absolute luminosity (~10^^ Lq) are required 
making the buried AGN scenario the most plausible choice (Soifer et al. 2003; their §6) 
In fact, known AGN populations (i.e. quasars) meet these requirements without difficulty, 
while for the starburst scenario, of order 100-1000 super star clusters must be concentrated 
in a very compact region and such a phenomenon has not been clearly confirmed from 



- 21 - 



observations. We thus believe that buried AGNs are the most natural explanation for the 
centrally-concentrated energy sources at ULIRGs' cores. The presence of lumi i ious b uried 
AGNs in many nearby ULIRGs agrees with the suggestion by iTakeuchi et al.l (120031 ) that 
the observed lump of the local infrared 60 fim luminosity function at the highest luminosity 
end is due to the AGN contribution. 



5.2.6. Combination of energy diagnostic methods 

Table 6 (column 7) summarizes the strengths of the detected buried AGN signatures in 
Spitzer IRS 5-35 fim spectra based on the three methods: (1) low PAH equivalent width; 
(2) large Tg j value; and (3) small r^g/rgy and t'{^/t'qj ratio. When buried AGN signatures 
in individual ULIRGs are consistently found using all or most of the methods, the ULIRGs 
are classified as very strong buried AGN candidates, marked with open double circles. When 
buried AGNs signatures are seen in the first method (particularly based on the observed 
EW6.2PAH value) and/or in the third method, then the ULIRGs are classified as strong AGN 
candidates (open circles). When the signatures are detected only in the second method 
and/or low observed EWu.spah or corrected EW6.2PAh values (first method), the ULIRGs 
are classified as possible buried AGN candidates (open triangles), because (1) a normal 
starburst nucleus obscured by foreground dust in the host galaxy (Figure Id) cannot be 
ruled out completely and (2) the 11.3 yum PAH emission strengths in normal starbursts 
have intrinsically larger scatter (§5.2.1). The fraction of ULIRGs with strong buried AGN 
signatures is 16/48 (33%; LINER + Hll-region), 10/28 (36%; LINER), and 6/20 (30%; HII- 
regions). When ULIRGs with possible buried AGN signatures are included, the fraction 
increases to 26/48 (54%; LINER + Hll-region), 19/28 (68%; LINER), and 7/20 (35%; HH- 
region). Our careful look at the infrared 5-35 /im spectra suggests that luminous buried 
AGNs reside in a much higher fraction of optical non-Seyfert ULIRGs than previous estimates 



flTaniguchi et al.lll999l : iLutz et al.lll999[ ) 



We stress that the above fraction of luminous buried AGNs is only for ULIRGs classified 
optically as LINERs and Hll-regions (i.e. non-Seyferts), which comprise 70% of ULIRGs. 
Given that 30% of ULIRGs show optical Seyfert signatures, indicative of luminous AGNs 
surrounded by dusty tori, and that the AGNs are e stimated to be energetically dominant in 



these optical Seyfert ULIRGs (IVeilleux et al.lll999bl ). the total fraction of ULIRGs possessing 



luminous AGNs is at least 50%, and possibly >65%, in the whole ULIRG sample at 2; < 
0.15, including both optical Seyfert and non-Seyfert objects. 

Finally, in all of the methods employed by us (§5.2.1-5.2.3 and 5.2.6), ULIRGs classified 
optically as LINERs tend to display signatures of luminous buried AGN more frequently than 
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those classified optically as Hll-regions. For starburst /buried AGN composite ULIRGs (Fig- 
ure Ic), it is likely that the optical LINER or Hll-region classifications are largely affected 
by the properties of the modestly-obscured starbursts at the exteriors of the ULIRG cores, 
rather than buried AGN-related emission, because optical observations can probe only the 
surfaces of dusty objects. If shocks or superwinds are important in the modestly-obscured 
starbursts, the optical classification will be LINERs. If emission from Hll-regions is dom- 
inant, the ULIRGs will be classed as Hll-regions. In this case, it is not obvious why two 
phenomena, a luminous buried AGN at the center and shocks in the surface starburst, are 
correlated. IVeilleux et al.l (Il995l . Il999al ) found that the emission probed in the optical is 
dustier in LINER ULIRGs than in Hll-region ULIRGs. In a dusty starburst, shock-related 
emission can be relatively important in the optical compared to the emission from Hll-regions 
themselves, resulting in optical LINER classification. When a luminous AGN is placed at 
the center of a less dusty starburst, classified optically as an Hll-region, the AGN emission 
is more easily detectable in the optical spectrum, making such an object an optical Seyfert. 
On the other hand, when a luminous AGN is placed at the center of a dusty starburst clas- 
sified optically as a LINER, the AGN emission is more elusive in the optical spectrum, so 
that such an object is classified as an optical non-Seyfert. The observed higher fraction of 
optically elusive buried AGN in optical LINER ULIRGs compared to Hll-region ULIRGs 
can qualitatively be explained by this scenario. 



5.2.7. Ah sorption- corrected luminosities of buried AGNs 

For ULIRGs which show buried AGN signatures and also small PAH equivalent widths, 
the observed fiuxes are largely ascribed to AGN-heated dust continuum emission. The result- 
ing luminosity is conserved at each temperature (Figure 2). We can quantitatively estimate 
the absorption-corrected intrinsic dust luminosity heated by the AGN (whole infrared), or 
the absorption-corrected AGN energetic radiation luminosity (X-ray - UV - optical), based 
on the observed fiuxes at Arest ~ 10 /xm and the dust column density toward the 10 /^m 
continuum emitting regions inferred from Tg y. Since in a buried AGN, the 10 fim continuum 
emission regions are dominated by fixed layers (§5.2.3; Figure 2), the foreground screen dust 
absorption model is applicable, and the dust extinction corrections are straightforward. 

To do this, we assume that Xg 7 and the extinction at Arest = 8 or 13 f im continuum just 



outsi de the 9.7 fim silicate feature (Acont) are related with Tg ^/Acont ~ 2.3 (IRieke &: Lebofsky 



19851 ). Table 7 summarizes the absorption-corrected AGN luminosities for selected ULIRGs 
with both strong buried AGN signatures and low PAH equivalent widths. The fiux attenu- 
ation of the 8 or 13 /im continuum outside the 9.7 /xm silicate feature ranges from a factor 
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of 1.2 (IRAS 01004-2237; r^^ ^ ~ 0.4) to 5 (IRAS 01298-0744; 7 ~ 4). The absorption- 
corrected AGN luminosities are 1-5.5 x 10^^ ergs s~^, which are of order ~1O^^L0 (= 3.8 x 
10^^ ergs s^^). 

The ambiguity of the derived absorption-corrected AGN luminosities come from (1) 
absolute flux calibration error of the Spitzer IRS spectra and slit loss of the compact AGN 
emission; (2) possibly different dust extinction curves from the Galactic one that we assume; 
and (3) a small contribution from PAH-emitting, modestly-obscured starbursts, which de- 
creases the Tg 7 values and increases the continuum flux level, compared to pure emission 
from the buried AGN. The first error is not significant, approximately 20% or so. The sec- 
ond uncertainty is more difficult to quantify, but our adoption of Tg y reduces the effects of 
dust grain size change (§5.2.2). For the third source of uncertainty, the effect of the decreased 
Tg 7 values is more important than that of the increased continuum flux level, because the 
AGN-produced flux at A 

rest ~ 9.7 yum is more highly attenuated. Assume for example that 
the modestly obscured starburst emission contributes 10% of the observed flux at Arest = 8 
and 13 /im, outside the 9.7 fim silicate feature. The buried AGN-produced continuum flux 
is then only 10% smaller than the observed flux. For ULIRGs with Tgy > 2 (i.e., buried 
AGN candidates we are now discussing), since the attenuation of the observed flux at Arest 
= 9.7 fim is a factor of >7 compared to the nearby continuum, more than 50% of the con- 
tinuum at Arest = 9.7 fim is produced by the starburst emission. The true, uncontaminated 
Tg y value toward the buried AGN is larger by a factor of >0.7 than the observed value, so 
that the correction factor of dust extinction for the AGN-heated dust continuum is larger by 
>30% (Tgy/Acont ~ 2.3). Therefore, the small starburst contamination generally results in 
the underestimate of the intrinsic absorption-corrected AGN luminosity. Despite these un- 
certainties, the rough agreement between the absorption-corrected buried AGN luminosities 
and the observed infrared luminosities of these ULIRGs (1-3 x makes it reasonable 

for us to argue that the putative AGNs could quantitatively account for the bulk of the 
observed infrared luminosities, at least for ULIRGs with strong buried AGN signatures in 
Spitzer IRS spectra. 



5.3. Comparison with infrared 3—4 jj,m spectra 



5.3.1. Ice absorption feature 



H2O ice absorption is also present at 3.1 /im (due to a stre tching mode) and has been 



detected in 3-4 /xm s pectra of many UL IRGs in this sample (jimanishi fc Maloneyl 12003 



Imanishi et al. 2006a; 



Risaliti et al.ll2006l ). Comparison of the 3.1 /im and 6.0 /im ice ab- 



sorption features can provide a better understanding of the properties of ice in ULIRGs. 
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Column 4 of Table 5 summarizes the detection or non-detection of the 3.1 fim ice 
feature. For the detections its optical depth is shown in column 5. Many ULIRGs have clear 
detections of both the 3.1 fim and 6.0 fim ice absorption features, reinforcing the picture 
that the energy sources in many ULIRGs are obscured by a large amount of ice-covered dust 
grains. 

Figure 12 compares the optical depths of the 3.1 /im (r^ i) and 6.0 /im (rg.o) H2O 
ice absorption features, for ULIRGs clearly showing both features. The tq q/t^ i ratios are 
several times larger than expected from laboratory data. Systematically larger tq q/t^^i ratios 
compared to the laborato ry prediction have previously been f ound in Galactic obscured stars, 
but the reason is unclear (iGibb et al.ll2000l : iKeane et al.ll200ll ). This trend may be common in 
astronomical obscured sources. Our poor understanding of ice absorption features hampers 
the search for buried AGNs with centrally-concentrated energy source geometry based on 
small T6,o/t3.i ratios, as was done based on the small t[^/t!^j ratio (§5.2.3). 



5.3.2. PAH emission features 



The 3.3 /im PAH emission feature assigned to the aromati c CH stretch ( A. 



lamandola et al. 



19891 ) was detected in the majority of this ULIRG sample (llmanishi et al.ll2006al ). Figure 
9 (d), (e) and (f) compare the rest-frame equivalent widths of the 3.3 /xm PAH feature 
(EW3.3PAH) with those of the 6.2 /im, 7.7 /im, and 11.3 /im PAH features. 

In Figure 9, the scatter is generally large in plots comparing PAH features at widely 
separated wavelengths. Contribution from PAH-free AGN continuum to an observed flux 
can vary to a larger degree with increasing wavelength separation, explaining the trend in 
Figure 9. 



5. 3. 3. Support for the strong dust temperature gradient 

In a buried AGN with a strong dust temperature gradient, dust extinction toward the 
3-4 /im continuum emission regions, estimated using 3-4 /im data, Av(3/xm), should be 
larger than that toward 10 /im continuum emission region estimated from the Tg ^ value 
(Figure 2). For ULIRGs with small contamination from PAH-emitting starbursts to their 
3-4 /im spectra (small EWs.spah values), Av(3/im) can be estimated with little ambiguity, 
from absorption optical depths at 3.4 /im (T34) by bare carbonaceous dust and at 3.1 /im 
(^3.1) by ice-covered dust deep inside molecular clouds. 

IRAS 08572+3915, 12127-1412, 17044+6720 show 3-4 /im spectra with very small PAH 
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contaminations. The t^a values are 0.8, 35, an d 0.15 for IRAS 08572+3915, 12127-1412 
17044+6720, respectively (llmanishi et al.ll2006a] ). Th e 4 value reflects the column density 



of bare carbonaceous dust (without an ice mantle) (jPendleton et al.l 1 19941 : llmanishi et al. 



19961 : iRawlings et al.ll2003l ). In the Galaxy, the t^j/t^^ ratio is surprisingly similar in dif- 



ferent directions {t^j/t3a = 14-17; Roche & Aitken 1984, 1985; Pendleton et al. 1994). 
Thus, a Tgj/r^A ratio substantially smaller than 14-17 can be a signature of a strong dust 
temperature gradient. Adopting the Xg 7 values shown in Table 4, we obtain Tg j/r^A = 5, 7, 
and 12 for IRAS 08572+3915, 12127-1412, 17044+6720, respectively. These three ULIRGs 
have previously been classifled as very strong buried AGN candidates, based on Spitzer IRS 
5-35 yum spectra (§5.2.6). The small Tg^/r^A ratios support the buried AGN hypothesis. 

For I RAS 08572+3915 and 17044+6720, no 3.1 fim absorption by ice-covered dust is 
detected (llmanishi et al.ll2006al ). suggesting that most of the dust toward the 3-4 fim con- 
tinuum emission sources is bare, without an ice mantle. However, for IRAS 12127—1412, the 
3.1 /im absorption feature is found with ra.i ~ 0.4, suggesting a large amount of ice-covered 
dust, in addition to bare dust grains. Since the 9.7 /im absorption is seen whether the silicate 
dust is bare or ice-covered, adding the column density of ice-covered dust to that of bare 
carbonaceous dust grains, will raise Av(3/im) even further, strengthening the argument for 
a strong dust temperature gradient in IRAS 12127—1412. 



5.3.4- Classification - AGN or starburst? 



Table 6 (column 8) summarizes the detection or non-detection, and the detection sig- 
niflcance in the case of detecti on, of buried AGN signatures based on ground-based 3-4 /im 
spectra (llmanishi et al.ll2006al ). The trend of buried AGN detection is generally consistent 
in individual ULIRGs, in that (1) very strong buried AGN candidates at 3-4 fim are also 
very strong buried AGN candidates in Spitzer IRS 5-35 /im spectra, and (2) ULIRGs with 
no AGN signatures at 3-4 /im usually show no AGN signatures at 5-35 /im. There are some 
ULIRGs which have detectable AGN signatures at 5-35 /im, but not at 3-4 /im, which are 
explained b y a larger contaminat i on from modes tly obscured starburst emission at a shorter 
wavelength (llmanishi et al.l |2004J : llmanishil |2006|) . 



Finally, llmanishi et al.l (l2006al ) have previously argued, based on 3-4 /im spectra, that 
the fraction of ULIRGs showing luminous buried AGN signatures is higher in LINER ULIRGs 
than Hll-region ULIRGs in this sample. This trend is exactly what we have seen in the 
Spitzer IRS 5-35 /im spectra discussed here (§5.2.1-5.2.6). 
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5.4. Buried AGNs and infrared colors 



Figure 13 compares the strength of the 9.7 fim sihcate absorption feature 



9.7) 



and 



infrared colors, Spitzer IRS 5.4 /xm to 35 /xm flux ratio (Figure 13a) and IRAS 25 yum to 60 
Hm flux ratio (Figure 13b). The infrared colors are independent of Tg 7 values, with no clear 
positive or negative correlation found. 

The IRAS 25 fim to 60 /xm flux ratio is often used to determine whether U LIRGs show 
cool (f25/f6o < 0.2) or warm (> 0.2) far-infrared colors (ISanders et al.l Il988bl ). Although 
many classical Seyfe rt-type AGNs surrounded by dusty tori show warm far-infrared colors 
(Ide Grijp et al.lll987l ). the bulk of ULIRGs possessing luminous buri ed AGN candid a tes (tq j 
> 2) belong to cool sources (Figure 13b). As already argued by llmanishi et al.l (l2006al ). 
it is quite reasonable that buried AGNs tend to show cooler infrared colors than optical 
Seyfert AGNs. A larger amount of dust around a central AGN in a buried AGN compared 
to a Seyfert AGN surrounded by a torus (§1) can naturally produce a cooler infrared color, 
because the contribution from an outer, cooler component becomes important, and makes 
the optical Seyfert signatures fainter or even vanish, because even the torus axis direction 
can be opaque to the bulk of the AGN's ionizing radiation. 



5.5. Future prospects 

Although infrared spectroscopic energy diagnostics enable us to discuss the overall prop- 
erties of ULIRGs statistically, follow-up observations of infrared-selected buried AGN can- 
didates are very important to cross-check the reliability of these methods. For example, 
we argue an exceptionally centrally-concentrated extreme starburst is unlikely, based on 
the emission surface brightness and absolute luminosity (§5.2.5). X-ray observations can be 
particularly powerful in distinguishing a buried AGN from an extreme starburst, because 
a luminous AGN is a much stronger X-ray emitter than any kind of starburst. Since the 
column densities in the cores of ULIRGs can be large enough to be Compton thick (X-ray ab- 
sorption Nh > 10^^ cm~^), high sensitivity X-ray observations at E > 10 keV are particularly 
desirable. 

In a strongly X-ray emitting buried AGN surrounded by dense gas and dust. X-ray 
dissociation regions (XDRs; Maloney et al. 1996) should develop, rather than the photo- 
dissociation regions (PDRs) usually seen in a starburst (a strong UV emitter) . Finding XDR 



signatures can thus be another powerful way to test the buried AGN scenario, llmanishi et al. 



(I2OO4J . l2006bl ) found that several infrared-selected buried AGN candidates tend to show 



millimeter molecular line flux ratios expected from XDRs. For one of the strongest infrared- 
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selected buried AGN candidates, IRAS 08572+3915, IShirahata et all (120071 ) detected CO 
(v=l-0) absorption features up to J = 17 in ground-based infrared 4.5-5 Min obse rvations. 
It is argued that such high-J lines can be explained only by XDRs mm, further 

supporting the buried AGN scenario. A further search for XDR signatures in ULIRGs may 
reinforce the buried AGN picture suggested from the infrared spectroscopy. 



6. Summary 

We presented Spitzer IRS 5-35 fim low-resolution (R ~ 100) spectra of ULIRGs aX z < 
0.15 in the IRAS 1 Jy sample. ULIRGs classified optically as non-Seyferts (LINERs and 
HII- regions), which comprise 70% of the 1 Jy sample at ^ < 0.15, were our prime targets to 
search for signatures of optically undetected AGNs deeply buried in dust along virtually all 
lines-of-sight, without well-developed narrow line regions. In total, 28 LINER and 20 HII- 
region ULIRGs (a complete sample) were observed. Three methods were mainly utilized: 
(1) the equivalent widths of the PAH emission features at 6.2 /xm, 7.7 /im, and 11.3 fim, (2) 
the absolute optical depth of the 9.7 fim silicate dust absorption feature, and (3) the 9.7 fim 
to 18 yum silicate dust absorption optical depth ratio. The first method traces the hardness 
of the radiation field (starburst-like soft or AGN-like hard). The latter two methods probe 
the geometry, i.e. whether the dust and energy sources are spatially well mixed (a normal 
starburst) or the energy source is more centrally concentrated than dust (a buried AGN). 
We found the following main conclusions. 

1. For the majority of ULIRGs, PAH emission features were detected, suggesting the 
presence of PAH-emitting normal starburst activity. However, the observed PAH to 
infrared luminosity ratios were systematically smaller than in modestly obscured (Ay 
< 20-30 mag) starburst galaxies. Taken at face value, the PAH-emitting starbursts 
with modest dust obscuration are energetically insignificant in ULIRGs. The dominant 
energy sources are (1) very highly obscured (Ay >> 30 mag) PAH-emitting starbursts 
with mixed dust/source geometry, whose PAH flux is attenuated, and/or (2) non- PAH- 
emitting AGNs. 

2. In 19 out of the 48 (40%) ULIRGs we studied, the observed PAH equivalent widths 
were found to be substantially lower than normal starburst galaxies at either of the 
6.2 /im, 7.7 /xm, or 11.3 /im PAH features. This fraction increased to 22/48 (46%) if 
we adopted the 6.2 /im PAH equivalent widths after correction for the ice and HAG 
absorption features. The highly obscured energy sources of these ULIRGs' cores were 
argued to be dominated by non-PAH-emitting AGNs, rather than highly-obscured 
PAH-emitting starbursts. 
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3. In 20/48 (42%) of ULIRGs, the observed optical depths of the 9.7 /im sihcate dust 
absorption feature (rg 7) were substantially larger than the upper limit determined by 
mixed dust /source geometry in normal starbursts. A non-PAH-emitting buried AGN 
with centrally-concentrated energy source geometry is the most natural explanation for 
the energy sources at the cores of these ULIRGs. Comparisons of the Xg 7 value with 
the PAH equivalent widths, and with the 11.3 /im to 6.2 fim PAH luminosity ratio, 
supported the buried AGN scenario in these ULIRGs. 

4. For selected ULIRGs with small PAH contamination (mostly ULIRGs with large Tg 7), 
the 18 fim to 9.7 yum silicate absorption optical depth ratios were often significantly 
smaller than the intrinsic ratio. A strong dust temperature gradient around a centrally- 
concentrated energy source geometry was strongly suggested, further supporting the 
buried AGN scenario. 

5. 16/48 (33%) of optical non-Seyfert ULIRGs showed strong signatures of luminous 
buried AGNs, and 26/48 (54%) ULIRGs displayed some signatures. Given that 30% of 
ULIRGs showed optical Seyfert signatures and were known to possess luminous AGNs 
surrounded by dusty tori, it was suggested that luminous AGNs are present in at least 
50%, and possibly >65%, of the complete ULIRG sample at z < 0.15. 

6. The fraction of ULIRGs showing luminous buried AGN signatures was always higher 
in optical LINER ULIRGs than Hll-region ULIRGs using any of the above methods. 
This was qualitatively explained by the scenario that the starburst regions probed in 
the optical are dustier in LINER ULIRGs than in Hll-region ULIRGs. 

7. For ULIRGs with strong buried AGN signatures and small PAH contamination, the 
absorption-corrected intrinsic luminosities of the buried AGNs were found to be quan- 
titatively comparable to IO^^Lq, accounting for the bulk of the observed infrared lu- 
minosities of ULIRGs. 

8. Most of the luminous buried AGN candidates were found in ULIRGs with cool far- 
infrared colors. 

9. The above conclusions drawn from the Spitzer IRS 5-35 /i m spectra were g enerall y 



consistent with those based on 3-4 jjm spectra obtained by llmanishi et al.l (j2006af ). 
Our main conclusions were that optically undetected luminous buried AGNs are surely 
present in a significant fraction of non-Seyfert ULIRGs. It is therefore very important to 
carefully and quantitatively investigate the energetic roles of optically elusive luminous 
buried AGNs, if we are to understand the true nature of the nearby ULIRG population. 
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A. The maximum Tgj value in a mixed dust/source geometry 

Figure la shows a schematic diagram of a typical mixed dust/source geometry. In this 
geometry, the observed flux I(r,^) is given by 

/(r,) = /o(A)xi^^, (Al) 

where /o(-^) is an unattenuated intrinsic flux and tx is the optical depth at each wavelength, 
which takes different values inside absorption features and nearby continua outside the fea- 
tures. 

We here consider radiation with an intrinsically smooth spectral energy distribution, 
namely Io(8 /xm) ~ Io(9.7 /im) ~ Io(13 /im), such as expected for dust emission in ULIRGs. 
The peak absorption optical depth of the 9.7 /im silicate feature, relative to the nearby 
continuum at 8 /zm and 13 /im (which corresponds to Tg 7 and not rg.y in our definition in 
§4.2), is 

Tcont 1 - e ^97.,! 

= In X , (A3) 

where tq^su and Tcont refer to the optical depths at the peak of the 9.7 /im silicate absorption 
feature and nearby continuum at 8 /im and 13 /im just outside the feature, respectively. The 
tqisU value is always larger than the Tcont value, because tqjsU is the superposition of the 
optical depth for continuum emission and the 9.7 fim silicate feature. 



In the case of the Galactic dust extinction curve derived by lRieke fc LebofskyI (119851 ). 



Tcont = 0.025/1.08 X Av (A4) 
T,Js^l = 0.087/1.08 X Av (A5) 

where, by definition. Ay = 1.08 x Ty. 

The absorption strength of the 9.7 /im silicat e feature, relative to Ay, is d ifferent by a 



factor of ~2 in different directions in the Galaxy (IRoche fc Aitkenlll984j . Il985l ). If we take 
this into account, 

Tcont = 0.025/1.08 X Ay (A6) 

TQ7sii = (0.05 ~ 0.1) X Av + 0.025/1.08 X Ay. (A7) 

The Tcont and Tgisu values can be arbitrarily large with increasing amount of dust (Ay) 
in the mixed dust /source geometry. The Tg 7 value is an increasing function of Ay, and the 
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first term in the equation A3 ( \_l-TZTii ) becomes unity for a very large Ay value. Thus, the 
upper limit to Tg y is 

T'g-j{upper — limit) = ln[^^^^^]. (A8) 

We obtain a stringent upper limit of T9 7 ~ 1.7 in this mixed dust/source geometry for any 
amount of dust. 

The Tcont/^v and Tgisii/A.Y ratios can vary substantially with varying dust grain size, 
because the wavelength difference between the optical V-band and infrared ~10 /xm is large. 
However, the upper limit of Tg ^ in the mixed dust/source geometry is dependent only on 
Tcont 1^97 sill which should be less sensitive to varying dust grain size (see Figure 6b), because 
their wavelengths are very close. Therefore, the above conclusion about the upper limit of 
Tg y in the mixed dust/source geometry should be robust to possible differences in dust grain 
size between the cores of ULIRGs and the Galactic diffuse interstellar medium. 
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Table 1. Observed ULIRGs at z < 0.15 and their IRAS-hased infrared emission properties 



Object 


Redshift 


fl2 




feo 


fioo 


log LiR 




Optical 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


L© 




Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(8) 


(9) 


TRAS 00188 — 0856 


0.128 


<0.12 


0.37 


2.59 


3.40 


12.3 


0.14 (C) 


LINER 


IRAS 00482—2721 


0.129 


<0.10 


<0.18 


1.13 


1.84 


12.0 


<0.16 (C) 


LINER 


IRAS 032504-lfiOfi 


0.129 


<0.10 


<0.15 


1.38 


1.77 


12.1 


<0.11 (C) 


LINER 


IRAS 04103—2838 


0.118 


0.08 


0.54 


1.82 


1.71 


12.2 


0.30 (W) 


LINER 


IRAS 08572+3915 


0.058 


0.32 


1.70 


7.43 


4.59 


12.1 


0.23 (W) 


LINER 


IRAS 09039+0503 


0.125 


0.07 


0.12 


1.48 


2.06 


12.1 


0.08 (C) 


LINER 


TRAS 0<3nfi+0SS4 


0.146 


<0.09 


<0.14 


1.09 


1.82 


12.1 


<'0 1 3 (C) 


LINER 


IRAS 09539+0857 


0.129 


<0.15 


<0.15 


1.44 


1.04 


12.0 


<0.11 (C) 


LINER 


IRAS 10378+1108 


0.136 


<0.11 


0.24 


2.28 


1.82 


12.3 


0.11 (C) 


LINER 


IRAS 10485—1447 


0.133 


<0.11 


0.25 


1.73 


1.66 


12.2 


0.14 (C) 


LINER 


IRAS 10494+4424 


0.092 


<0.12 


0.16 


3.53 


5.41 


12.1 


0.05 (C) 


LINER 


IRAS 11095-0238 


0.106 


0.06 


0.42 


3.25 


2.53 


12.2 


0.13 (C) 


LINER 


IRAS 11130-2659 


0.136 


<0.09 


0.20 


1.21 


1.24 


12.1 


0.17 (C) 


LINER 


IRAS 12112+0305 


0.073 


0.12 


0.51 


8.50 


9.98 


12.3 


0.06 (C) 


LINER 


IRAS 12127—1412 


0.133 


<0.13 


0.24 


1.54 


1.13 


12.1 


0.16 (C) 


LINER 


IRAS 12359—0725 


0.138 


0.09 


0.15 


1.33 


1.12 


12.1 


0.11 (C) 


LINER 


IRAS 13335—2612 


0.125 


<0.13 


<0.14 


1.40 


2.10 


12.1 


<0.10 (C) 


LINER 


IRAS 14252—1550 


0.149 


<0.09 


<0.23 


1.15 


1.86 


12.2 


<0.20 (C) 


LINER 


IRAS 14348-1447 


0.083 


0.07 


0.49 


6.87 


7.07 


12.3 


0.07 (C) 


LINER 


IRAS 15327+2340 (Arp 220) 


0.018 


0.48 


7.92 


103.33 


112.40 


12.1 


0.08 (C) 


LINER 


IRAS 16090-0139 


0.134 


0.09 


0.26 


3.61 


4.87 


12.5 


0.07 (C) 


LINER 


IRAS 16468+5200 


0.150 


<0.06 


0.10 


1.01 


1.04 


12.1 


0.10 (C) 


LINER 


IRAS 16487+5447 


0.104 


<0.07 


0.20 


2.88 


3.07 


12.1 


0.07 (C) 


LINER 


IRAS 17028+5817 


0.106 


<0.06 


0.10 


2.43 


3.91 


12.1 


0.04 (C) 


LINER 


IRAS 17044+6720 


0.135 


<0.07 


0.36 


1.28 


0.98 


12.1 


0.28 (W) 


LINER 


IRAS 21329-2346 


0.125 


0.05 


0.12 


1.65 


2.22 


12.1 


0.07 (C) 


LINER 


IRAS 23234+0946 


0.128 


<0.06 


0.08 


1.56 


2.11 


12.1 


0.05 (C) 


LINER 


IRAS 23327+2913 


0.107 


<0.06 


0.22 


2.10 


2.81 


12.1 


0.10 (C) 


LINER 


IRAS 00091-0738 


0.118 


<0.07 


0.22 


2.63 


2.52 


12.2 


0.08 (C) 


HII 


IRAS 00456-2904 


0.110 


<0.08 


0.14 


2.60 


3.38 


12.2 


0.05 (C) 


HII 


IRAS 01004-2237 


0.118 


0.11 


0.66 


2.29 


1.79 


12.3 


0.29 (W) 


HII 


IRAS 01166-0844 


0.118 


0.07 


0.17 


1.74 


1.42 


12.1 


0.10 (C) 


HII 


IRAS 01298-0744 


0.136 


<0.12 


0.19 


2.47 


2.08 


12.3 


0.08 (C) 


HII 


IRAS 01569-2939 


0.141 


<0.11 


0.14 


1.73 


1.51 


12.2 


0.08 (C) 


HII 


IRAS 02411+0353 


0.144 


<0.08 


0.22 


1.37 


1.95 


12.2 


0.16 (C) 


HII 


IRAS 10190+1322 


0.077 


<0.07 


0.38 


3.33 


5.57 


12.0 


0.11 (C) 


HII 


IRAS 11387+4116 


0.149 


0.12 


<0.14 


1.02 


1.51 


12.2 


<0.14 (C) 


HII 


IRAS 11506+1331 


0.127 


<0.10 


0.19 


2.58 


3.32 


12.3 


0.07 (C) 


HII 


IRAS 13509+0442 


0.136 


0.10 


<0.23 


1.56 


2.53 


12.2 


<0.15 (C) 


HII 


IRAS 13539+2920 


0.108 


<0.09 


0.12 


1.83 


2.73 


12.0 


0.07 (C) 


HII 


IRAS 14060+2919 


0.117 


<0.10 


0.14 


1.61 


2.42 


12.1 


0.09 (C) 


HII 


IRAS 15206+3342 


0.125 


0.08 


0.35 


1.77 


1.89 


12.2 


0.20 (C) 


HII 


IRAS 15225+2350 


0.139 


<0.07 


0.18 


1.30 


1.48 


12.1 


0.14 (C) 


HII 


IRAS 16474+3430 


0.111 


<0.13 


0.20 


2.27 


2.88 


12.1 


0.09 (C) 


HII 
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Table 1 — Continued 



Object 




Redshift 


fl2 


f25 


feo 


fioo 


log LiR 




Optical 








(Jy) 


(Jy) 


(Jy) 


(Jy) 


L© 




Class 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


IRAS 20414- 


-1651 


0.086 


<0.65 


0.35 


4.36 


5.25 


12.2 


0.08 (C) 


HII 


IRAS 21208- 


-0519 


0.130 


<0.09 


<0.15 


1.17 


1.66 


12.0 


<0.13 (C) 


HII 


IRAS 22206- 


-2715 


0.132 


<0.10 


<0.16 


1.75 


2.33 


12.1 


<0.09 (C) 


HII 


IRAS 22491- 


-1808 


0.076 


0.05 


0.55 


5.44 


4.45 


12.1 


0.10 (C) 


HII 



Note. — Col.(l): Object name. Col.{2): Redshift. Col.(3)-(6): fi2. f25. ko, and fipp are IRAS 
fluxes at 12 ^m, 25 fim, 60 fim, and 100 ^m, respectively, taken from [Kim &: Sanders! lll998l ). Col. (7): 
Decimal logarithm of infrared (8—1000 fj,m) luminosity in units of solar luminosity (L p)), calculated with 
LiR = 2.1 X 1039 X D{Mpc)2 X (13.48 x /12 + 5.16 x /25 + 2.58x /go + Zioo) ergss"! ISand crs & Mirabel 
Il996r) . Since the calculation is based on our adopted cosmology, the infrared luminosities slightly (<10%) 
differ from the values shown in Kim & Sanders (1998, their Table 1, column 15). For sources that have 
upper limits in some IRAS bands, we can derive upper and lower limits of the infrared luminosity, by 
assuming that the actual flux is the IRAS-wpper limit and zero value, respectively. The difference of the 
upper and lower values is usually very small, less than 0.2 dex. We assume that the infrared luminosity 
is the average of these values. Col. (8): IRAS 25 ^m to 60 fim flux ratio. ULIRGs with f25/f60 < 0-2 
and > 0.2 are classified as cool and warm sources (denoted as "C " and " W"), respectively llSanders et al.l 
Il988d) . Col. (9): Optical spectral classification bv lveilleux et al.l lll999ah . 
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Table 2. Spitzer IRS observing log 



Object 


PID 


Date 




Integration time 


[sec] 








[UT] 


SL2 


SLl 


LL2 


LLl 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IRAS 00188-0856 


105 


2003 Dec 17 


240 


240 


180 


180 


IRAS 00482-2721 


3187 + 2306 


2005 Jul 7 + 2004 Jul 18 


240 


240 


240 


240 


IRAS 03250+1606 


3187 + 2306 


2005 Feb 11 


240 


240 


240 


240 


IRAS 04103-2838 


3187 + 2306 


2005 Feb 10 + 2004 Aug 10 


240 


240 


240 


240 


IRAS 08572+3915 


105 


2004 Apr 15 


84 


84 


84 


84 


IRAS 09039+0503 


3187 + 2306 


2005 Apr 18 + 2004 Nov 17 


240 


240 


240 


240 


IRAS 09116+0334 


2306 


2005 Apr 21 


240 


240 


240 


240 


IRAS 09539+0857 


3187 + 2306 


2005 Jun 5 + 2005 May 31 


240 


240 


240 


240 


IRAS 10378+1108 


105 


2005 Jun 8 


240 


240 


180 


180 


IRAS 10485-1447 


3187 + 2306 


2005 May 23 + 2005 Jan 4 


240 


240 


480 ^ 


480 a 


IRAS 10494+4424 


2306 


2004 Nov 17 


240 


240 


240 


240 


IRAS 11095-0238 


105 


2005 Jun 7 


240 


120 


240 


120 


IRAS 11130-2659 


2306 


2005 Jul 12 


240 


240 


240 


240 


IRAS 12112+0305 


105 


2004 Jan 4 


84 


84 


120 


120 


IRAS 12127-1412 


3187 + 2306 


2005 Jun 30 


240 


240 


240 


240 


IRAS 12359-0725 


2306 


2005 Jun 30 


240 


240 


240 


240 


IRAS 13335-2612 


3187 + 2306 


2005 Feb 15 + 2004 Jul 17 


240 


240 


240 


240 


IRAS 14252-1550 


2306 


2004 Jul 17 


240 


240 


240 


240 


IRAS 14348-1447 


105 


2004 Feb 7 


120 


120 


120 


120 


Arp 220 


105 


2004 Feb 29 


84 


84 


60 


60 


IRAS 16090-0139 


105 


2005 Aug 6 


120 


120 


180 


180 


IRAS 16468+5200 


2306 


2004 Jul 14 


480 


480 


480 


480 ^ 


IRAS 16487+5447 


2306 


2004 Jul 17 


240 


240 


240 


240 


IRAS 17028+5817 


2306 


2004 Jul 17 


240 


240 


240 


240 


IRAS 17044+6720 


2306 


2004 Jul 17 


240 


240 


240 


240 


IRAS 21329-2346 


3187 + 2306 


2004 Nov 16 + 2004 Oct 23 


240 


240 


240 


240 


IRAS 23234+0946 


3187 + 2306 


2004 Dec 13 + 2004 Dec 8 


240 


240 


240 


240 


IRAS 23327+2913 


2306 


2004 Dec 8 


240 


240 


240 


240 


IRAS 00091-0738 


3187 + 2306 


2005 Jun 30 + 2004 Dec 8 


240 


240 


240 


240 


IRAS 00456-2904 


3187 + 2306 


2005 Jul 14 + 2004 Jul 17 


240 


240 


240 


240 


IRAS 01004-2237 


105 


2004 Jan 4 


120 


120 


120 


120 


IRAS 01166-0844 


3187 + 2306 


2005 Jan 3 + 2005 Jan 15 + 2004 Jul 17 


480 


480 


480 ^ 


480 


IRAS 01298-0744 


105 


2005 Jul 14 


240 


240 


120 


120 


IRAS 01569-2939 


2306 


2004 Jul 18 


240 


240 


240 


240 


IRAS 02411+0353 


2306 


2005 Jan 14 


240 


240 


240 


240 


IRAS 10190+1322 


3187 + 2306 


2005 May 22 + 2004 Dec 12 


240 + 240 240 + 240 


480 * 


480 * 


IRAS 11387+4116 


2306 


2005 Jan 11 


240 


240 


240 


240 


IRAS 11506+1331 


3187 + 2306 


2005 May 25 + 2005 Jan 3 


240 


240 


240 


240 


IRAS 13509+0442 


2306 


2004 Jul 17 


240 


240 


240 


240 


IRAS 13539+2920 


2306 


2005 Feb 7 


240 


240 


240 


240 


IRAS 14060+2919 


2306 


2004 Jul 16 


240 


240 


240 


240 


IRAS 15206+3342 


105 


2004 Jun 24 


120 


120 


180 


180 


IRAS 15225+2350 


2306 


2005 Feb 7 


240 


240 


240 


240 


IRAS 16474+3430 


2306 


2004 Jul 14 


240 


240 


240 


240 
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Table 2 — Continued 



Object 




PID 


Date 


Integration time 


[sec] 








[UT] 


SL2 


SLl 


LL2 


LLl 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IRAS 20414- 


-1651 


105 


2004 May 14 


120 


120 


120 


120 


IRAS 21208- 


-0519 


3187 + 2306 


2004 Nov 13 + 2004 Oct 22 


240 <= 


240 ° 


480 


a 4gQ a 


IRAS 22206- 


-2715 


3187 + 2306 


2004 Nov 15 + 2004 Oct 25 


240 


240 


480 


a a 


IRAS 22491- 


-1808 


105 


2004 Jun 24 


120 


120 


120 


120 



^Although separate spectroscopy of double nuclei (|Kim et al.ll20o3 ) was proposed, the same coor- 
dinate was observed twice, probably because IRS blue peak up camera (13-18.5 iim) pointed to the 
same position. 

We combined SL spectra of double nuclei with ~3" separation l|Kim et al.ll2002h . because both 
nuclei are very faint. The combined SLl spectrum is smoothly connected to the LL2 spectrum. 

'^We combined SL spectra of double nuclei with ^^5" separation l|Kim et al.ll2002h . because the 
fainter northern nucleus is very faint. The combined SLl spectrum is smoothly connected to the LL2 
spectrum. 

SL spectra of the eastern and western nuclei with ~4" separation llKim et al.ll2002h are analyzed 
separately, because both nuclei have similar flux levels. The combined SLl spectrum is smoothly 
connected to the LL2 spectrum. 

®Only the brighter northern nucleus llKim et al.ll200 j) was analyzed, because (1) the fainter south- 
ern nucleus with ~7" separation is very faint, (2) it may not be fully covered with the LL spectrum, 
and (3) the SLl spectrum of the northern nucleus is smoothly connected to the LL2 spectrum. 

Note. — Col.{l): Object name. Col.{2): PID number: 105 (PI = J. Houck), 2306 (PI = M. 
Imanishi), and 3187 (PI = S. Veilleux). For sources denoted with "3187 -I- 2306", SL and LL 
spectra were taken by observations with PID = 3187 and 2306, respectively. Col. (3): Observing 
date in UT. Col. (4): Net on-source integration time for SL2 spectroscopy in sec. Col. (5): Net on- 
source integration time for SLl spectroscopy in sec. Col. (6): Net on-source integration time for LL2 
spectroscopy in sec. Col. (7): Net on-source integration time for LLl spectroscopy in sec. 



Table 3. PAH emission 



Object 


EW6.2PAH 


EW7.7PAH 


EW11.3PAH 


EW3.3PAH 


L6.2PAH 


L7.7PAH 


L11.3PAH 


L6.2PAh/LiR 


Lii.3pah/Lir 




[nm] 


[nm] 


[nm] 


[nm] 


10*2 [ergs s-i] 


10"'^ [ergs s^i] 


lO-'^ [ergs s-i] 


[X 10-3] 


[x 10-3] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IRAS 00188-0856 


85 (~40) 


305 


290 


50 


3.5 


17 


3.4 


0.4 


0.4 


IRAS 00482-2721 


260 


710 


465 




1.8 


4.8 


0.9 


0.5 


0.2 


IRAS 03250+1606 


325 


740 


570 


80 


7.0 


18 


4.6 


1.6 


1.1 


IRAS 04103-2838 


165 


365 


120 




8.4 


21 


6.6 


1.4 


1.1 


IRAS 08572+3915 


<15 


<100 


<30 


<5 


<1.9 


<18 


<0.5 


<0.4 


<0.1 


IRAS 09039+0503 


325 (~200) 


625 


580 


95 


3.7 


11 


2.2 


0.8 


0.5 


IRAS 09116+0334 


360 


730 


630 


75 


9.4 


27 


7.4 


1.9 


1.5 


IRAS 09539+0857 


145 


405 


390 


65 


2.5 


12 


1.3 


0.6 


0.3 


IRAS 10378+1108 


110 (~50) 


450 


210 


40 


2.2 


14 


1.5 


0.3 


0.2 


IRAS 10485-1447 


225 (~150) 


400 


150 


80 


3.7 


11 


1.1 


0.6 


0.2 


IRAS 10494+4424 


345 (~250) 


715 


495 


110 


6.1 


19 


3.1 


1.1 


0.6 


IRAS 11095-0238 


20 {~10) 


250 


110 


150 


1.0 


17 


0.9 


0.2 


0.1 


IRAS 11130-2659 


95 (~60) 


490 


280 




1.9 


13 


1.0 


0.4 


0.2 1 


IRAS 12112+0305 


390 


565 


425 


100+160 


6.1 


15 


3.2 


0.8 


0.4 


IRAS 12127-1412 


5 (~3) 


95 


20 





0.7 


11 


0.6 


0.1 


0.1 ° 


IRAS 12359-0725 


220 


445 


165 


75 


3.1 


8.8 


1.1 


0.6 


0.2 


IRAS 13335-2612 


310 


700 


610 




6.6 


17 


5.1 


1.5 


1.2 


IRAS 14252-1550 


225 


685 


320 


60+60 


3.6 


12 


2.3 


0.6 


0.4 


IRAS 14348-1447 


260 (~110) 


540 


500 


70+110 


5.9 


19 


3.6 


0.8 


0.5 


Arp 220 


215 (~120) 


410 


315 


80 


1.5 


5.1 


0.8 


0.3 


0.2 


IRAS 16090-0139 


110 (~70) 


265 


325 


75 


8.1 


30 


5.6 


0.7 


0.5 


IRAS 16468+5200 


230 


470 


435 


120 


2.9 


11 


2.0 


0.7 


0.5 


IRAS 16487+5447 


250 


580 


570 


70 


2.8 


9.6 


2.2 


0.5 


0.4 


IRAS 17028+5817 


375 


800 


595 


120 


4.6 


14 


2.8 


0.9 


0.6 


IRAS 17044+6720 


55 


165 


30 


10 


3.9 


13 


1.0 


0.8 


0.2 


IRAS 21329-2346 


230 


635 


370 


50 


3.0 


12 


1.3 


0.6 


0.3 


IRAS 23234+0946 


220 


695 


445 


75 


3.1 


11 


2.0 


0.7 


0.4 


IRAS 23327+2913 


255 


485 


175 


45 


2.5 


6.9 


1.4 


0.6 


0.3 


IRAS 00091-0738 


50 (~30) 


350 


255 




1.3 


15 


1.5 


0.2 


0.3 


IRAS 00456-2904 


385 


750 


545 




11 


27 


7.6 


1.9 


1.4 


IRAS 01004-2237 


45 


190 


8 




3.1 


19 


0.5 


0.4 


0.1 


IRAS 01166-0844 


30 


430 


190 




0.5 


8.5 


0.7 


0.1 


0.2 


IRAS 01298-0744 


10 (-6) 


380 


155 




0.7 


37 


1.1 


0.1 


0.1 



Table 3 — Continued 



\jD]ecz 


i!j W6.2PAH 


JiW7.7PAH 


T7\A/, . ...^ . 


VV3.3PAH 


L6.2PAH 


a 

L7.7PAH 


L11.3PAH 


L6.2PAh/LiR 


Lii.3PAh/Lir 




[nm] 


[nm] 


[nm] 


[nm] 








[v in— 3i 


[X iU J 


('-) 




{•^) 


(A\ 

I*; 






(1) 








IRAS 01569-2939 


115 


555 


255 




3.1 


20 


2.0 


0.5 


0.3 


IRAS 02411+0353 


300 


695 


375 




13 


38 


7.9 


2.1 


1.3 


IRAS 10190+1322E 


375 


680 


510 


95 


4.4 


12 


3.0 


1.2 


0.8 


IRAS 10190+1322W 


375 


745 


685 


50 


4.6 


11 


3.6 


1.2 


0.9 


IRAS 11387+4116 


265 


680 


475 


75 


4.2 


15 


3.9 


0.7 


0.7 


IRAS 11506+1331 


195 (~120) 


540 


305 


95 


8.5 


26 


3.4 


1.1 


0.4 


IRAS 13509+0442 


490 


795 


590 


135 


9.7 


22 


4.7 


1.5 


0.7 


IRAS 13539+2920 


355 


755 


600 


85 


7.6 


22 


4.4 


1.9 


1.1 


IRAS 14060+2919 


340 


860 


530 


150 


8.6 


27 


6.5 


2.0 


1.5 


IRAS 15206+3342 


225 


525 


230 


55 


12 


32 


7.7 


2.0 


1.3 


IRAS 15225+2350 


145 (~100) 


290 


130 


40 


5.0 


16 


2.0 


1.0 


0.4 


IRAS 16474+3430 


390 


820 


575 


105 


9.2 


27 


5.8 


1.7 


1.1 


IRAS 20414-1651 


330 


705 


465 


75 


3.2 


10 


1.7 


0.5 


0.3 


IRAS 21208-0519 


505 


810 


420 


100 


6.0 


14 


2.7 


1.5 


0.7 *^ 


IRAS 22206-2715 


210 


975 


565 




3.1 


13 


2.7 


0.6 


0.5 


IRAS 22491-1808 


325 


695 


460 




3.9 


12 


2.6 


0.8 


0.5 



^We regard flux excess at Arost = 7.3—8.1 /im above an adopted continuum level as 7.7 /im PAH emission, to reduce the effects of the strong 9.7 fim silicate dust absorption 
feature. Our definition is different from those in many previous literatures. 



Note. — Col.(l): Object name. Col. (2): Rest-frame equivalent width of the 6.2 /ira PAH emission. The values in parentheses are approximate ones after correction for the ice 
and HAG absorption features. Col. (3): Rest-frame equivalent width of the 7.7 fim PAH em ission. Col. (4): Rest-frame equivalent width of the 11.3 /im PAH emission. Col. (5): 
Rest-frame equivalent width of the 3.3 ^m PAH emission, taken from llmanishi et al.l ll2006al ). — : no data. Col. (6): Luminosity of the 6.2 fim PAH emission in units of 10*^ ergs 
. Col. (7): Luminosity of the 7.7 /im PAH emission in units of lO*^ ergs s~^. Col. (8): Luminosity of the 11.3 /^m PAH emission in units of lO''^ ergs s '^ Col.(9): The 6. 2 
/xm PAH to infrared luminosity ratio in units of 10^"^. The ratio for normal starbursts with modest dust obscuration (Ay < 20—30 mag) is ~3.4 X 10~^ jPeeters et al.ll2004l ). 
Col. (10): The 11.3 /tm PAH to infrared luminosity ratio in units of 10~^. The ratio for normal starbursts with modest dust obscuration (Av < 20-30 mag) is ~1.4 X 10"'^ 
l lSoifer et al.ll2002h . 
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Table 4. Optical depth of the 9.7 /xm and 18 /xm silicate dust absorption features 



Object 


119.7 


/ 

^9.7 


^18 


^11 


as/ ^9. 7 


// / f 

^18/^9.7 


(1) 


(2) 


(3) 


/ ,1 \ 

(4) 


(5) 


(6) 


(7) 


IRAS 00188—0856 


3.1 


2.5 


0.55 


0.65 


0.22 


0.26 


IRAS 00482—2721 


2.6 


2.1 


0.6 


0.75 


0.29 


0.36 


IRAS 03250+1606 


1.8 


1.6 










IRAS 04103—2838 


0.5 


0.6 










IRAS 08572+3915 


4.1 


3.8 


0.9 


1.0 


0.24 


0.26 


IRAS 09039+0503 


2.6 


2.2 


0.9 


1.0 


0.41 


0.45 


IRAS 09116+0334 


1.5 


1.1 










IRAS 09539+0857 


3.8 


3.5 


1.3 


1.6 


0.37 


0.46 


IRAS 10378+1108 


2.6 


2.4 


0.45 


0.55 


0.19 


0.23 


IRAS 10485—1447 


3.0 


3.0 


0.8 


0.95 


0.27 


0.32 


IRAS 10494+4424 


2.2 


1.7 










IRAS 11095-0238 


3.5 


3.4 


0.8 


0.9 


0.23 


0.26 


IRAS 11130-2659 


3.5 


3.1 


1.0 


1.2 


0.32 


0.39 


IRAS 12112+0305 


2.3 


1.7 










IRAS 12127-1412 


3.0 


2.5 


0.6 


0.7 


0.24 


0.28 


IRAS 12359-0725 


2.2 


1.8 


0.5 


0.55 


0.28 


0.31 


IRAS 13335-2612 


1.6 


1.4 










IRAS 14252-1550 


1.6 


1.4 










IRAS 14348-1447 


2.8 


2.0 










Arp 220 


3.7 


3.2 


0.95 


1.1 


0.30 


0.34 


IRAS 16090-0139 


3.0 


2.6 


0.7 


0.8 


0.27 


0.31 


IRAS 16468+5200 


3.0 


2.5 


0.65 


0.8 


0.26 


0.32 


IRAS 16487+5447 


2.4 


1.8 










IRAS 17028+5817 


2.0 


1.5 










IRAS 17044+6720 


1.7 


1.8 


0.4 


0.5 


0.22 


0.28 


IRAS 21329-2346 


2.7 


2.2 










IRAS 23234+0946 


2.2 


1.8 










IRAS 23327+2913 


1.5 


1.3 










TR nOOQI —ClJ'^R 
irxrto uuuy± — u/oo 


3.2 


3.2 


1.05 


1.25 


0.33 


0.39 


IRAS 00456-2904 


1.5 


1.2 










IRAS 01004-2237 


0.4 


0.4 










IRAS 01166-0844 


3.2 


3.0 


0.75 


0.9 


0.25 


0.30 


IRAS 01298-0744 


4.3 


4.0 


1.1 


1.35 


0.27 


0.34 


IRAS 01569-2939 


3.2 


2.8 


1.0 


1.2 


0.36 


0.43 


IRAS 02411+0353 


1.3 


1.0 










IRAS 10190+1322 


1.7 


1.3 










IRAS 11387+4116 


1.6 


1.1 










IRAS 11506+1331 


2.6 


2.3 


0.7 


0.85 


0.30 


0.37 


IRAS 13509+0442 


2.0 


1.6 










IRAS 13539+2920 


1.8 


1.6 










IRAS 14060+2919 


1.5 


1.0 










IRAS 15206+3342 


0.9 


0.6 










IRAS 15225+2350 


2.5 


2.2 


0.65 


0.75 


0.29 


0.34 


IRAS 16474+3430 


1.8 


1.5 










IRAS 20414-1651 


2.4 


1.9 
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Table 4 — Continued 



Object 

(1) 


(2) 


(3) 


(4) 


(5) 


'''is/'''').? '''is/'''').? 
(6) (7) 


IRAS 21208-0519 


1.5 


1.4 








IRAS 22206-2715 


2.4 


1.7 








IRAS 22491-1808 


2.0 


1.5 









'^To estimate the optical depths of silicate dust absorption features, SL 
spectra of the eastern and western nuclei are summed. 

Note. — Col.(l): Object name. Col. (2); Tg.Y is an optical depth of 
the 9.7 //m silicate dust absorption feature, against a power-law continuum 
determined from data points at Arest = 5.6 //m, 7.1 /im and Aobs = 34-35 
/um, shown as dashed lines in Figures 3 and 4. Col. (3): Tg y is an optical 
depth of the 9.7 ^m silicate dust absorption feature, against a power-law 
continuum determined from data points at Arost = 7.1 fim and 14.2 fim, 
shown as dotted lines in Figures 3 and 4. Once the continuum levels are 
fixed, the uncertainty of Tg y is <5% for ULIRGs with large Tg ^ values (>2) 
and can be ~10% for ULIRGs with small Tg y. Col. (4): T(g is an optical 
depth of the 18 /im silicate dust absorption feature, against a power-law 
continuum determined from data points at Arcst = 14.2 fim and 24 /im, 
shown as dotted lines in Figures 3 and 4. Once the continuum is fixed, the 
statistic uncertainty of t-[j^ is <10%, because the value is estimated only for 
ULIRGs with clearly detectable 18 ixm silicate absorption features. Col. (5); 
T(g is an optical depth of the 18 /im silicate dust absorption feature, against 
a power-law continuum determined from data points at Arest = 14.2 ;um and 
29 ^m, shown as dashed-dotted lines in Figures 3 and 4. Once the continuum 
is fixed, the statistic uncertainty is <10%. Col. (6): t-J^/t/j ^ ratio, only for 
ULIRGs with clearly detectable 18 fiiu silicate absorption. The uncertainty is 
<10%, once the continuum is fixed. Col. (7): T{g/Tg y ratio, only for ULIRGs 
with clearly detectable 18 /xm silicate absorption. 
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J. Ct ULXD O . 


Ice absorpt 


1 on 




Ohiprt 


6.0 nm 


T6.0 


3.1 jum 


'^3.1 




(2) 


(3) 






IRAS 00188-0856 


o 


0.9 


o 


1.8 


IRAS 00482-2721 


X 








IRAS 03250+1606 


X 




o 


0.8 


IRAS 04103-2838 


X 








IRAS 08572+3915 


X 




X 




IRAS 09039+0503 


O 


0.6 


o 


0.8 


IRAS 09116+0334 


X 




o 


0.4 


IRAS 09539+0857 


X 




X 




IRAS 10378+1108 


O 


0.8 


o 


0.6 


IRAS 10485-1447 


O 


0.5 


o 


0.8 


IRAS 10494+4424 


O 


0.4 


o 


1.0 


IRAS 11095-0238 


o 


0.5 


X 




IRAS 11130-2659 


o 


0.5 






IRAS 12112+0305 


A 




X 




IRAS 12127-1412 


o 


0.4 


O 


0.4 


IRAS 12359-0725 


X 




o 


0.5 


IRAS 13335-2612 


X 








IRAS 14252-1550 


X 




o 


0.7 


IRAS 14348-1447 


O 


0.9 


o 


0.5 


Arp 220 


O 


0.6 






IRAS 16090-0139 


o 


0.5 


o 


0.8 


IRAS 16468+5200 


A 




X 




IRAS 16487+5447 


A 




o 


1.2 


IRAS 17028+5817 


A 




o 


0.6 


IRAS 17044+6720 


X 




X 




IRAS 21329-2346 


A 




o 


1.0 


IRAS 23234+0946 


X 




X 




IRAS 23327+2913 


X 




X 




IRAS 00091-0738 


O 


0.5 






IRAS 00456-2904 


X 








IRAS 01004-2237 


X 








IRAS 01166-0844 


X 








IRAS 01298-0744 


O 


0.5 






IRAS 01569-2939 


X 








IRAS 02411+0353 


X 








IRAS 10190+1322E 


X 




X 




IRAS 10190+1322W 


X 




X 




IRAS 11387+4116 


X 




X 




IRAS 11506+1331 


O 


0.5 


O 




IRAS 13509+0442 


X 




X 




IRAS 13539+2920 


X 




X 




IRAS 14060+2919 


X 




X 




IRAS 15206+3342 


X 




X 




IRAS 15225+2350 


O 


0.5 


O 


0.4 


IRAS 16474+3430 


X 




O 


0.8 
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Table 5 — Continued 



Object 




6.0 iim 


f"6.0 


3.1 fim 




(1) 




(2) 


(3) 


(4) 


(5) 


IRAS 20414- 


-1651 


X 




X 




IRAS 21208- 


-0519 


X 




X 




IRAS 22206- 


-2715 


X 








IRAS 22491- 


-1808 


X 









Note. — Col.(l): Object name. Col. (2): Detection or 
non-detection of the 6.0 /im ice absorption feature. Q- 
clearly detected. A: possibly detected. X: non-detected. 
The non-detection includes two cases; (i) the absorption 
feature is indeed weak or absent in high signal-to-noise ra- 
tios spectra, or (ii) the presence of ice absorption is not 
clear due to relatively large scatter in spectra or to the dif- 
ficulty in separating the ice absorption and nearby PAH 
emission features. Col. (3): Optical depth of the 6.0 fim ice 
absorption feature. — : not measurable. Since data scatter 
in spectra is not the same for individual ULIRGs, there is 
no fixed detection threshold. The uncertainty is ~0.1 for 
detected sources. Col. (4): Detection or non-detection of 
th e 3.1 fim ice absorption feature in the 2.8-4.1 fim spectra 
of llmanishi et al] l l2006al ). Q- clearly detected. X: non- 
detected. — : no data. Col. (5): Optical depth of the 3.1 
fim ice absorption feature. — : not measurable because (i) 
it is weak, (ii) data scatter is large in spectra, and (ill) 2.8- 
4.1 fim spectra are lacking. Again there is no fixed detection 
threshold. 
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Table 6. Buried AGN signatures 



Object 






Spitzer IRS 5- 


-35 iMa 






3-4 nm 




E^Vg 2PAH 


EW7_7PAH 


-^^^ll.SPAH 


t' 
9.7 


T-gradient 


Total 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


IRAS 00188-0856 


O 


X 


X 


o 


o 


® 


O 


IRAS 00482-2721 


X 


X 


X 


o 


X 


A 




IRAS 03250+1606 


X 


X 


X 


X 


X 


X 


O 


IRAS 04103-2838 


O 


X 


o 


X 


X 


O 




IRAS 08572+3915 


O 


o 


O 


o 


O 


® 


® 


IRAS 09039+0503 


X 


X 


X 


o 


X 


A 


o 


IRAS 09116+0334 


X 


X 


X 


X 


X 


X 


o 


IRAS 09539+0857 


O 


X 


X 


o 


X 


o 


X 


IRAS 10378+1108 


O 


X 


X 


o 


O 


o 


o 


IRAS 10485-1447 


x(0) 


X 


O 


o 


A 


A 


o 


IRAS 10494+4424 


X 


X 


X 


X 


X 


X 


o 


IRAS 11095-0238 


O 


X 


O 


o 


O 


@ 


X 


IRAS 11130-2659 


o 


X 


X 


o 


X 


O 


X 


IRAS 12112+0305 


X 


X 


X 


X 


X 


X 


X 


IRAS 12127-1412 


o 


o 


o 


o 


O 


® 


® 


IRAS 12359-0725 


X 


X 


O 


X 


A 


A 


o 


IRAS 13335-2612 


X 


X 


X 


X 


X 


X 




IRAS 14252-1550 


X 


X 


X 


X 


X 


X 


o 


IRAS 14348-1447 


x(0) 


X 


X 


X 


X 


A 


o 


Arp 220 


x(0) 


X 


X 


O 


X 


A 


X 


IRAS 16090-0139 


O 


X 


X 


o 


A 


O 


o 


IRAS 16468+5200 


X 


X 


X 


o 


A 


A 


X 


IRAS 16487+5447 


X 


X 


X 


X 


X 


X 


o 


IRAS 17028+5817 


X 


X 


X 


X 


X 


X 


o 


IRAS 17044+6720 


O 


o 


o 


X (-) 


O 


@ 


@ 


IRAS 21329-2346 


X 


X 


X 


O 


X 


A 


o 


IRAS 23234+0946 


X 


X 


X 


X 


X 


X 


X 


IRAS 23327+2913 


X 


X 


O 


X 


X 


A 


X 


IRAS 00091-0738 


O 


X 


X 


O 


X 


O 




IRAS 00456-2904 


X 


X 


X 


X 


X 


X 




IRAS 01004-2237 


O 


o 


O 


X(-)- 


X 


® 




IRAS 01166-0844 


O 


X 


O 


O 


O 


@ 




IRAS 01298-0744 


O 


X 


O 


O 


X 


O 




IRAS 01569-2939 


o 


X 


X 


O 


X 


O 




IRAS 02411+0353 


X 


X 


X 


X 


X 


X 




IRAS 10190+1322 


Xb 


Xb 


Xb 


X 


X 


X 


X 


IRAS 11387+4116 


X 


X 


X 


X 


X 


X 


X 


IRAS 11506+1331 


x(0) 


X 


X 


O 


X 


A 


O 


IRAS 13509+0442 


X 


X 


X 


X 


X 


X 


X 


IRAS 13539+2920 


X 


X 


X 


X 


X 


X 


X 


IRAS 14060+2919 


X 


X 


X 


X 


X 


X 


X 


IRAS 15206+3342 


X 


X 


X 


X 


X 


X 


X 


IRAS 15225+2350 


O 


X 


o 


O 


X 


o 


O 


IRAS 16474+3430 


X 


X 


X 


X 


X 


X 


O 
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Table 6 — Continued 



Object 








Spitzer IRS 5-35 fj.m 






3-4 fim 






EW6.2PAH 


EW7.7PAH 


EWii,3PAH 


T-' 


T-gradient 


Total 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


IRAS 20414- 


-1651 


X 


X 


X 


X 


X 


X 


X 


IRAS 21208- 


-0519 


X 


X 


X 


X 


X 


X 


X 


IRAS 22206- 


-2715 


X 


X 


X 


X 


X 


X 




IRAS 22491- 


-1808 


X 


X 


X 


X 


X 


X 





''The small Tg ^ value is due to weak obscuration of a luminous AGN, rather than the normal starburst 
nature with mixed dust /source geometry. 

''PAH equivalent widths are derived for both nuclei separately. No AGN signatures are found in both nuclei. 

Note. — Col.{l): Object name. Col. (2): Buried AGN signatures based on the low equivalent width of the 6.2 
/im PAH emission (EW6.2PAH < 180 nm) (§5.2.1). Q- present. X: none. When the corrected EW6.2PAH value 
becomes less than 180 nm, the symbol "O" is shown in parentheses. Col. (3): Buried AGN signatures based on 
the low equivalent width of the 7.7 /im PAH emission {EW7,7pah < 230 nm) (§5.2.1). Q- present. X: none. 
Col. (4): Buried AGN signatures based on the low equivalent width of the 11.3 fj,m PAH emission (EWn 3pah 
< 200 nm) (§5.2.1). Q: present. X: none. Col. (5): Buried AGN signatures based on the large Tg y value (>2) 
(§5.2.2). O- present. X: none. Col. (6): Buried AGN signatures based on the small T^g/Tg y and ''""g/rg ^ 
ratios. O- present. A: possibly present. X: none. Considering the uncertainty of the ratio with <10% (Table 
4), ULIRGs showing Tjg/rg ^ < 0.27 and T"g/rg < 0.32 are taken as sources with detectable signatures of 
strong dust temperature gradients (§5.2.3). Col. (7): Buried AGN signatures from combined methods in Col. 
(2)-(6) based on Spitzer IKS spectra. @: very strong. Q: strong. A: possibly present. X: none. Please see 
§5.2.6 for more details. Col. (7): Buried AGN signatures based on infrared 3—4 fim spectroscopic diagnostic by 
llmanishi et al. very strong. Q- strong. X: none. — ; no data. 
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Table 7. Absorption-corrected AGN luminosity 



Object 


L(AGN) [10^5 g^gs g-i] 


(1) 


(2) 


IRAS 00188-0856 


1.5 


IRAb 04103—2838 


1.5 


IRAS 08572+3915 


5.5 


IRAS 10378+1108 


1.0 


IRAS 10485-1447 


1.0 


IRAS 11095-0238 


2.5 


IRAS 11130-2659 


1.0 


IRAS 12127-1412 


2.5 


IRAS 16090-0139 


2.0 


IRAS 17044+6720 


2.0 


IRAS 00091-0738 


1.5 


IRAS 01004-2237 


1.5 


IRAS 01166-0844 


1.0 


IRAS 01298-0744 


3.5 


IRAS 15225+2350 


1.5 



Note. — Col.(l): Object name. Col.(2): 
Absorption-corrected intrinsic luminosity of the 
AGN energetic radiation (X-ray - UV - optical) 
in units of 10^^ ergs s~^. 
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PAH + continuum 





PAH-free continuum 



Buried AGN 



PAH + continuum 
(weal<ly obscured) 



dust 
in a host galaxy) 



PAH-free continuum 
(highly obscured) 



Buried AGN 




PAH + continuum 
(absorbed) 



4^ 





















PAH + continuum 
(weakly obscured) 



PAH-free continuum 
(highly obscured) 



Extreme starburst 
(Hll-regions only) 

Fig. 1. — (a): Geometry of energy sources and dust in a normal starburst. The open star 
symbols indicate "stars" in a starburst. The energy sources (stars) and dust are spatially 
well mixed. Large equivalent width PAH emission is observed, (b): Geometry of the energy 
source and dust in a buried AGN. The energy source (a mass-accreting supermassivc black- 
hole) is more centrally concentrated than the surrounding dust. No PAH emission is found. 
(c): A buried AGN/starburst composite. The starburst surrounds the central buried AGN. 
The observed spectrum is a superposition of PAH + continuum emission from the starburst 
and PAH-free continuum from the AGN. Since the buried AGN is more highly obscured than 
the surrounding starbursts, the starburst emission generally makes a strong contribution to 
an observed flux even if the starburst is energetically insignificant, (d): A normal starburst 
nucleus (mixed dust/source geometry) obscured by a large amount of foreground dust in a 
host galaxy. This geometry can produce a strong dust absorption feature whose optical depth 
is larger than the maximum threshold obtained in the mixed dust /source geometry, (e): An 
exceptionally centrally-concentrated extreme starburst, whose emitting volume is predomi- 
nantly occupied by Hll-regions. Such an extreme starburst produces no PAH emission, like 
a buried AGN. 
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Fig. 2. — Schematic diagram of a strong dust temperature gradient around a centrally- 
concentrated energy source (the central small filled circle). Inner dust close to the central 
energy source has higher temperature than dust far away from the energy source. Assuming 
approximately blackbody radiation, dust at lOOOK dust, close to the innermost dust subli- 
mation radius, shows an emission peak at Arest ~ 3 fim. More distant 300K dust and even 
150K dust have emission peaks at Arest ~ 10 fim and ~20 fim, respectively. This is a very 
simplified picture, and in an actual case, dust with other temperatures can also contribute 
to the observed fiux at each wavelength. For any reasonable dust radial density distribution, 
it is always true that the contribution from inner, hotter dust emission is larger at a shorter 
wavelength in this wavelength range. When we look at the energy source from the right 
side, an observed fiux is dominated by dust at the foreground side (right), rather than dust 
with the same temperature at the distant side (left), due to fiux attenuation, particularly for 
outer cooler dust. Thus, the relation Av(3yum) > Av(lOyum) > Av(20/im) should hold. Dust 
extinction toward the 3 /im continuum emitting region, Av(3/im), can be estimated from 
the optical depths of absorptioi i features at 3.4 p,m due to bare carbonac eous dust and at 
3.1 fim due to ice-covered dust (jimanishi et al.ll200ll . l2006al : Ilmanishill2006l ). The extinction 
toward the 10 /im and 20 fim continuum emitting regions, Av(10/im) and Av(20/im), can be 
estimated from the optical depths of 9.7 fim and 18 fim silicate dust absorption features in 
Spitzer IRS spectra, respectively. Although a spherical geometry is illustrated here for sim- 
plicity, the essence is unchanged as long as a centrally-concentrated compact energy source 
is obscured by dust along our line-of-sight. 
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Fig. 3. — Spectra of ULIRGs classified optically as LINERs. The abscissa and ordinate 
are, respectively, the observed wavelength in /im and fiux F^, in Jy, both shown as decimal 
logarithmic scale. For all objects, the ratio of the uppermost to lowermost scale in the 
ordinate is a factor of 1000, to show the variation of the overall spectral energy distribution. 
Dashed line: power-law continuum determined from data points at Arcst = 5.6 /xm, 7.1 fim, 
and Aobs = 34-35 /im. Dotted line: power-law continuum determined from data points at 
Kcst = 7.1 fim and 14.2 fira for 9.7 fim silicate absorption, and at Aj-cst = 14.2 fira and 24 fim 
for 18 iJ,m silicate absorption. Dashed-dotted line: power-law continuum determined from 
data points at Arest = 14.2 /im and 29 /im, used only for ULIRGs showing strong 18 /xm 
silicate dust absorption features. 



-54- 




10 20 
wavelength (pm) 



10 20 
wavelength (jj,m) 



10 20 
wavelength (jj,m) 



-55- 




5 10 20 5 10 20 ^5 10 20 

wavelength (pm) wavelength (jim) wavelength (/im) 




^5 10 20 5 10 20 5 10 20 



wavelength (/tim) wavelength (jj,m) wavelength (jj,m) 




5 10 20 5 10 2( 

wavelength (pm) wavelength (jj,m) 



Fig. 4. — Spectra of ULIRGs classified optically as Hll-regions, shown in the same way as 
Figure 3. Symbols are also the same as Figure 3. For IRAS 10190+1322, SL spectra of the 
eastern and western nuclei are combined. 
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Fig. 5. — Spectra of all ULIRGs at Aobs = 5.2-14.5 /xm, to zoom in on the PAH emission 
features and display them in detail. The abscissa and ordinate are, respectively, observed 
wavelength in //m and flux in Jy, both shown at hnear scale. For IRAS 10190+1322, the 
spectra of the eastern and western nuclei are shown separately. The expected wavelengths 
of the 6.2 /xm, 7.7 //m, and 11.3 /xm PAH emission features, as well as 6.0 /xm H2O ice 
absorption feature, are indicated as down arrows. The mark "(?)" is added, when detection 
is unclear. The solid lines are adopted continuum levels to estimate the strengths of the 
detected PAH emission features. 



10 20 30 40 

wavelength (/Um) 



10 20 30 40 

wavelength (/um) 



Fig. 6. — Silicate dust absorption profile, (a): Solid lines are the observed silica te ab- 
sorption profiles of the Galactic stars GCS 3-1 and WR 98a (jChiar fc Tielend l2006l ). The 
dashed line is the theoretical ab sorption profile of "astronomical silicate" grains with 0.1 
/xm in size ( jPraine fc Led Il984j ). All profiles are normalized at the 9.7 /xm absorption 
peak. The silicate absorption profile extends from 8 / zm to >30 /iin. f b): Theoret- 
ical absorption profile of "astronomical silicate" grains f Praine Sz Leel 1984 ) of different 
sizes, 1 fim (dotted line), 0.1 /im (dashed line), and 0.01 /xm (solid line). Adopted from 



http://www.astro.princeton.edu/~draine/dust/dust.diel.html The ordinate is in arbitrary 



units, but is proportional to the actual optical depth. In this figure, only absorption is taken 
into account, and scattering is neglected, because scattering is negligible in this wavelength 
range (see the above webpage). 
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Fig. 7. — IRS spectra at Aobs — 5.2-9 //m for ULIRGs displaying clear 6.0 iim H2O absorption 
feature (marked with "6.0/xm ice absorption" in the first 15 plots). Spectra of the next 
five ULIRGs [marked with "6.0/im ice absorption (?)"] show possible signatures of the 
6.0 fim ice absorption feature. Spectra of the last two ULIRGs, IRAS 17044+6720 and 
00456—2904, marked with "no 6.0;um ice absorption" are shown as examples with undetected 
ice absorption features. The abscissa is observed wavelength in /xm in linear scale, and the 
ordinate is flux in Jy in decimal logarithmic scale. 
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Fig. 8. — (a): Comparison of the rest-frame equivalent width of the 6.2 /xm PAH emission 
(abscissa) and 6.2 /xm PAH to infrared luminosity ratio (ordinate). Filled triangle: ULIRGs 
classified optically as LINERs. Open circle: ULIRGs classified optically as Hll-regions. The 
horizontal dotted line indicates the ratio for normal starburst galaxies with modest dust 
obscuration (L6.2pah/Lir = 3.4 x 10~^; see §5.1). In the upper-left small area separated by 
the dashed fine, there are no data points, (b): The abscissa is the rest- frame equivalent width 
of the 11.3 //m PAH emission and the ordinate is the 11.3 //m PAH to infrared luminosity 
ratio. The horizontal dotted line indicates the ratio for normal starburst galaxies with modest 
dust obscuration (Lii.3pah/Lir = 1-4 x 10~^; see §5.1). The larger scatter of Lh.spah/Lir 
in this Figure than for L6.2pah/Lir in Figure 8a may be due to intrinsically larger scatter of 
the 11.3 /im PAH strengths in normal starburst galaxies (see §5.2.1). 



-64- 




400 

'6.2PAH [nrn] 



200 

EW, 



400 
6.2PAH [""^] 



600 




E 



a. 

CM 

(d 



_ g 



O 

o 

U3 



O 
O 



O 

o 

CM 



(d) 
- ▲ 


LINER 1 o o 


1 ' ' 




o 


HII 1 








CAT A ^ 


o 






A aO A 

O A 






▲ 










1 1 I 1 1 1 1 1 1 I 1 







50 



100 



150 



200 



7.7PAH 



EW 



3.3PAH 



[nm] 



o 
o 
o 



-(e) A LINER I 
O HII I 



-I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I- 



o 



' A^A 



CP A o 



O 



O' 



_i I I I I I I I I I I I I I I I I i_ 



o 
o 
o 



o 
o 

? " 

o 

i " 

o 
o 



o 
o 

CM 



50 



EW 



100 

3.3PAH 



150 



200 



—I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 

(f) ^ LINER I 
O HII I 



o 



aaOao O 

O 

A o 



[nm] 



O 



50 



100 
'3.3PAH 



150 



200 



EW, ,n*u [nm] 



Fig. 9. — (a): Comparison of tlic rest-frame equivalent widths of PAH emission features in nm. 
The abscissa is the 6.2 jim PAH and the ordinate is the 7.7 iim PAH emission. Fihed triangle: 
ULIRGs classified optically as LINERs. Open circle: ULIRGs classified optically as HH-regions. 
In the upper-left small area separated with the dashed line, there are no data points, (b): The 
abscissa is the 6.2 /im PAH and the ordinate is the 11.3 /im PAH emission, (c): The abscissa is 
the 7.7 ^m PAH and the ordinate is the 11.3 /um PAH emission, (d): The abscissa is the 3.3 nm 
PAH and the ordinate is the 6.2 ^um PAH emission, (e): The abscissa is the 3.3 /im PAH and the 
ordinate is the 7.7 //m PAH emission, (f): The abscissa is the 3.3 //m PAH and the ordinate is 
the 11.3 //m PAH emission. In the last three plots, the average EW3.3PAH value of both nuclei is 
plotted for IRAS 12112+0305 and 14348-1447 (Table 3). 
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Fig. 10. — (a): Comparison of the strength of the 9.7 //m sihcate dust absorption (rg y) (ab- 
scissa) and the rest-frame equivalent widths of the 6.2 //m PAH emission feature (ordinate). 
Filled triangle: ULIRGs classified optically as LINERs. Open circle: ULIRGs classified op- 
tically as Hll-regions. In the upper-right small area separated with the dashed line, there 
are no data points, (b): The ordinate is the rest-frame equivalent width of the 11.3 fim PAH 
emission feature. In both plots, the two outliers IRAS 01004—2237 and 17044+6720 (weakly 
obscured AGNs; §5.2.4) are indicated as "IR01004" and "IR17044", respectively. For IRAS 
10190-1-1322, the average value of both nuclei (Table 3) are plotted for PAH equivalent widths 
in the ordinate. 
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Fig. 11. — Comparison of the strength of the 9.7 fim sihcate dust absorption (rg^) (ab- 
scissa) and 11.3 /xm to 6.2 //m PAH luminosity ratio (ordinate). Filled triangle: ULIRGs 
classified optically as LINERs. Open circle: ULIRGs classified optically as Hll-regions. In 
the upper-left small area separated with the dashed line, there are no data points. For IRAS 
10190+1322, the 11.3 fim to 6.2 fim PAH luminosity ratio in the ordinate is derived by 
combining spectra of both nuclei. IRAS 08572+3915 is not plotted, because both the 6.2 
fim and 11.3 fim PAH emission luminosities are upper limits (Table 3). 
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Fig. 12. — Comparison of the optical depths of the H2O ice absorption features at 3.1 /im 
(abscissa) and 6.0 /im (o rdinate). The das hed hne is the expected ratio from a laboratory 
H2O ice sample at 10 K (iKeane et al.ll200ll ). 
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Fig. 13. — Comparison of the strength of the 9.7 /xm sihcate dust absorption and infrared 
continuum colors, (a): The abscissa is the Tq j value and the ordinate is the 5.4 /xm to 35 /im 
flux ratio measured from the Spitzer IRS spectra in Figures 3 and 4. Filled triangle: ULIRGs 
classified optically as LINERs. Open circle: ULIRGs classified optically as Hll-regions. In 
the upper-right small area separated with the dashed line, there are no data points, (b): 
The ordinate is the IRAS 25 to 60 fim flux ratio, taken from T able 1. The horizonta l 
dotted line indicates the separation between warm and cool ULIRGs (ISanders et al.lll988bl ). 



